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1 Introduction 

1.1 Unique observable 

Analysis of the azimuthal anisotropy resulting from non-central nuclear collisions 
appears to be one of the most informative directions in studying the nature and 
properties of matter created in high energy nuclear collisions. Anisotropics in par- 
ticle momentum distributions relative to the reaction plane, often referred to as 
anisotropic collective flow or event anisotropics, have been in use for a few decades, 
starting from the first Berkeley Bevalac experiments. Azimuthal anisotropics at- 
tracted even more attention when the so called in-plane elliptic flow, first suggested 
as a signature of collective flow in relativistic nuclear collisions by OUitrault HI, 
was experimentally observed at the Brookhaven Alternate Gradient Synchrotron 
(AGS) IHIS, and later at the CERN Super Proton Synchrotron (SPS) [4|. At the 
Brookhaven Relativistic Heavy Ion Collider (RHIC) the observation of large ellip- 
tic flow f?\ is considered one of the most important discoveries which lead to the 
concept of the strongly interacting Quark Gluon Plasma (sQGP). Anisotropic flow 
will be among the first results at the CERN Large Hadron Collider (LHC) heavy-ion 
program. 

The main interest in anisotropic flow is due to its sensitivity to the system proper- 
ties very early in its evolution. The origin of anisotropics in the particle momentum 
distributions lies in the initial asymmetries in the geometry of the system. Because 
the spatial asymmetries rapidly decrease with time, anisotropic flow can develop 
only in the first fm/c. Based on this, one can conclude that anisotropic flow must 
be sensitive to the particle interactions very early in the system evolution, informa- 
tion usually available only via weakly interacting probes. In this sense, anisotropic 
flow is a unique hadronic observable providing direct information about the stage 
where the QGP may be the main player. Constituent rescattering is by far the most 
common explanation of anisotropic flow. Although possibilities of a different origin 
of elUptic flow have been discussed (e.g. partonic structure of the nuclei 161, color 
dipole orientation fTl, or direct anisotropy in particle emission from the Color Glass 
Condensate (CGC) ISj^) we do not consider them here. 

Previous review papers ifTOl [TTI on collective flow in heavy-ion collisions have 
presented results from accelerators at lower energy than RHIC. In this review we 
concentrate on results obtained in recent years, with the data coming mostly from 
RHIC. The field of anisotropic flow is growing with new data and theoretical results 
appearing rapidly. In this review we concentrate on the general results and inter- 
pretation, the status of the field, and major unresolved issues. The idea is that an 
interested reader, and we expect many graduate students and young researches to 
be among them, could not only appreciate the achievements of this field, but also 
identify interesting problems and be ready to start working in those directions. For 
that reason we also give a rather detailed presentation of the "technical" side of flow 
measurements, discussing advantages and disadvantages of different methods and 
associated systematic uncertainties in the results. Though for the real details we re- 
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fer to the original papers, we hope that the information presented here could provide 
a good basis to get involved. One will find that, unfortunately, the systematic uncer- 
tainties in flow measurements are still rather large, up to 10-15%, and often more. 
We identify two directions for future flow measurements, one being large statistics 
to try to better understand systematics, and the other being measurement of flow of 
rare particles. 



1.2 Definitions: flow and nonflow, the reaction and participant 
planes. 

The reaction plane is spanned by the vector of the impact parameter and the beam 
direction. Its azimuth is given by f Rp. The particle azimuthal distribution measured 
with respect to the reaction plane is not isotropic; so it is customary to expand it in 
a Fourier series [12 J: 

where the v„ = (cos[n(0,- — fRp)]) coefficients are used for a quantitative charac- 
terization of the event anisotropy, and the angle brackets mean an average over all 
particles in all events. The sine terms are not present because of symmetry with re- 
spect to the reaction plane, vi is referred to as directed flow, and V2 as elliptic flow 
(see Fig. [TJ. Radial flow in this paper refers to radial in the transverse plane. The 
Vn coefficients are functions of rapidity and transverse momentum, and as such they 
are often referred to as n''' harmonic differential flow. By integrated flow we mean 
the values of the v„ coefficients averaged over transverse momentum and rapidity. 




Fig. 2 The definitions of the Reaction Plane 
Fig. 1 Diagrams of elliptic and directed flow. and Participant Plane coordinate systems. 
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The reaction plane angle can not be directly measured in high energy nuclear 
collisions, but can be estimated from the particle azimuthal distribution event-by- 
event. Then the different harmonic flow coefficients are reconstructed from two or 
many particle azimuthal correlations. This introduces uncertainty in the analysis, 
discussed in more detail in the methods section, as the azimuthal correlations are not 
determined solely by anisotropic flow but have other contributions, usually referred 
to as nonflow and in case of two-particle correlations quantified by parameters 5„: 

(cosK0,-0;)]) = (v'') + 5„. (2) 

Anisotropic flow can fluctuate event to event, both in magnitude and direction 
even at fixed impact parameter We desciihe flow fluctuations by 

(^l={^i)-{^'nf- (3) 

One of the important sources of flow fluctuations are fluctuations in the initial geom- 
etry of the overlapping region due to the random nature of the interaction between 
constituents of the two nuclei. The participants are those constituents which partake 
in the primary interaction. The principal axis of the participant zone can deviate 
from the reaction plane. Fig.[2]shows the axes in the participant coordinate system, 
compared to the reaction plane system. It is important to distinguish between flow 
values measured in these two systems; the values in the reaction plane system be- 
ing always smaller than in the participant plane system: V2,pp > t^'2.RP- We discuss 
flow fluctuations due to fluctuations in the initial participant zone geometry in more 
detail in section U.3.41 



2 Experimental methods 
2.1 Event plane method 

In the standard event plane method ||3l [T5]| one estimates the azimuthal angle of the 
reaction plane from the observed event plane angle determined from the anisotropic 
flow itself. This is done for each harmonic, «, of the Fourier expansion. The event 
flow vector Qn is a 2d vector in the transverse plane: 

Qn.x = ^w,cos(«0,) = Q„cos(n'f^„), 

Qn.y = X^vv,sin(«0,) = Q„sin(«f„), (4) 

i 

where the sum goes over all particles / used in the event plane calculation. The 
quantities 0/ and w, are the lab azimuthal angle and weight for particle /, where 
for odd harmonics w,(— y) — ~Wi{y). The optimal choice for is to approximate 
Vn{pT,y)- Since often v„{pj',y) almost linearly increases with pT, the transverse 
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momentum is a common choice as a weight. The event plane angle is the azimuthal 
angle of Q„ calculated as 

f„ = arctan2(e„,y,e„,.0/n, (5) 

where arctan2 is a C language mathematical function. 

The observed v„ is the n''^ harmonic of the azimuthal distribution of particles with 
respect to this event plane: 

^°n\PT,y) - (cos[«(0,-'f^,)]), (6) 

where angle brackets denote an average over all particles in all events with their az- 
imuthal angle 0, in a given rapidity and pj momentum space bin at a fixed centrality. 
To remove auto-correlations one has to subtract the Q-vector of the particle of in- 
terest from the total event Q-vector, obtaining a to correlate with the particle. To 
avoid binning problems one should store the cosine directly in a profile histogram, 
rather than making a histogram of — f^, and then obtaining the mean cos. 

Since finite multiplicity limits the estimation of the angle of the reaction plane, 
the v„ have to be corrected for the event plane resolution for each harmonic given 
by 

^„ =(cos[«(f„-'f^Rp)]), (7) 

where angle brackets denote an average over a large event sample. The final flow 
coefficients are 

,,obs 

This equation should be applied in a narrow centrality bin. For a wide centrality bin, 
one should average the results from the narrow bins weighted with the multiplicity 
of the bin, since v„ is a particle-wise average. 

The reaction plane resolution depends on the multiplicity of particles used to 
define the flow vector and the average flow of these particles via the resolution 
parameter |fT3l[T4l[T5l: 



Z = v„ VM (9) 

= Zexp(-;tV2) (/(,-i)/2UV2) + Vi)/2(zV2)), (10) 
where / is the modified Bessel function. (Note that the definition of parameter % in 



Ref. ifTSl was larger by V2.) The dependence of Eq. (10 1 on is shown for the case 
of = 1 in Fig.|3] To estimate the event plane resolution one divides the full event up 
into two independent sub-events fTS' TT^ of equal multiplicity. Since the sub-events 
are positively correlated because each is correlated with the reaction plane, the event 
plane resolution for the sub-events is just the square-root of this correlation: 



(11) 
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where A and B denote the two subgroups of particles. Given ^„,sub, the solution for 



% in Eq. (10 1 is done by iteration. The full event plane resolution is obtained using 



Eq. ( 10 1 from the resolution of the sub-events by 



(12) 



because yM and the full event has twice as many particles as the sub-events. 
In the low resolution (< 0.5) linear region of the graph ioxk— 1, ^fuii ~ V^^sub- 




Fig. 3 The event plane resolution as a function of v„,\/M. The harmonic number of the correlation 
n is an integer k times the harmonic number m of the event plane. 



There may be reasons not to use the Q-vector of the full event, but to correlate 
particles from one sub-event with the Q-vector of the other sub-event. In this case 
the resolution of the sub-event plane should be used, and the particle of interest is 
automatically not included in the Q-vector Criteria which have been used for divid- 
ing the event into sub-events are: random, pseudorapidity, charge, and combinations 
of these. Using sub-events separated in pseudorapidity is a good way of reducing 
contributions from short-range correlations, as flow is a large scale effect. 

To remove acceptance correlations from an imperfect detector, one must first 
make the Q-vector in Eqs. (j4]) isotropic in the laboratory, both for the sub-events 
and the full event (if needed). Three methods have been used ifTSll for this flattening 
of the event plane azimuthal distribution: 

1 . Phi Weighting - one weights each particle with the inverse of the azimuthal dis- 
tribution of the particles averaged over many events. 

2. Recentering - one subtracts from the Q-vector of each event, the Q-vector aver- 
aged over many events. 

3. Shifting - one fits the non-flat distribution of averaged over many events with 
a Fourier expansion and calculates the shifts for each event f^, necessary to force 
a flat distribution on average. 
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By "many events" we mean a large enough sample to obtain good averages, but 
small enough to avoid shifts in the beam position and/or detector response as a func- 
tion of time. The first method is more intuitive, while the second is more practical 
because it guarantees zero average Q-vector It is also less sensitive to strong varia- 
tions in acceptance and can deal with "holes" in the detector. If either of the first two 
methods is not sufficient, then the third method may also be used. However, it should 
be pointed out that only the n''^ harmonic of the flattened distribution needs to be 
small when = 1 ; that is when one is not dealing with mixed harmonics. Calculating 
the distribution of — "f^, and dividing by that for mixed events has no advantage 
over the phi-weight method. A complete rigorous treatment of acceptance effects 
can be achieved in the cumulant approach using generating functions 1 18 1 and sim- 
ilarly, but in a somewhat more transparent way, in the scalar product method 1 19 1. 

An event plane determined from harmonic m allows one to study the flow of 
harmonics n = km, where k is an integer. In Eqs. Q and (|5]l for the event plane, 
n is the harmonic number of the event plane, but the n in Eq. (j6|l is the harmonic 
number of the correlation, which must be an integral multiple k of the event plane 
harmonic. The case of A; > 1 is called the mixed harmonics method. It was used 
widely at the AGS and SPS, where in the fixed target setting the detectors usually 
cover well the region of rapidity where directed flow is large. At RHIC it is mostly 
used to study higher (n > 4) harmonics relative to elliptic flow, e.g. V4{EP2}. This is 
useful because elliptic flow at RHIC is very strong near midrapidity. However, it is 
also true that determining the event plane from directed flow of neutrons in a Zero 
Degree Calorimeter, allows one to greatly suppress nonflow effects in elliptic flow 
measurements at midrapidity. 

The resolution from Eq. ( 10 1 is lower in the mixed harmoruc method (see Fig. |3]l 
for the case A; = 2, because of the added difficulty of resolving higher harmonics. 
However, the advantages are that nonflow correlations are greatly suppressed be- 
cause one is correlating two different harmonics of the collective flow, which is 
essentially a three-particle correlation. Also, one can determine the sign of the cor- 
relation harmonic relative to the event plane harmonic. The sign of vi is defined to 
be positive for nucleons at large positive rapidity. 

Examples of results from the event plane method are shown in Fig. [4] vi is an 
odd function of rapidity, with the proton flow in the opposite direction to the pion 
flow. V2 is an even function of rapidity, peaking at midrapidity. As a function of pj 
all curves go to zero at zero pj. 

Unless specified otherwise, all flow values presented in this paper will be from 
the full event plane method with the event plane determined near mid-rapidity. 



2.2 Two and many particle correlations 



The pair-wise correlation method ll20ll . is based on the fit of the two-particle az- 
imuthal distribution to that expected from anisotropic flow: 
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— — c<(l+^2y2cosM(/.)) (13) 

where all pairs of particles in a given momentum region are correlated. No event 
plane is used. Acceptance correlations are removed to first order by dividing by the 
mixed event distribution. The harmonic coefficients are small because they are the 
squares of the flow coefficients. This equation is only for integrated quantities, but 
normally the integrated values are obtained by averaging the differential quantities. 

The two-particle cumulant method differs from the previous one only in that 
instead of the fit to the two-particle distribution, it calculates the coefficients directly 
as 

y„{2}2 = (cos[«(0i - 02)]) = (m«jm:,2> (14) 

for all pairs of particles, where m„ = e'"'^ is a particle's unit flow vector. The differ- 
ential flow in the scalar product method ||2TI differs from the event plane method of 
Eq. ([8]l by using the magnitude of the flow vector as a weight 



where m„ ,■ is the unit vector of the /'''-particle (which is not included in Q„) and a 
and b are two subevents. The resulting statistical errors are slightly smaller than with 
the standard event plane method. If Q„ is replaced by its unit vector, this reduces to 
the standard method. The differential flow from the two-particle cumulant method 
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is the same as Eq. [TSjwith each Q-vector weighted by one over its muhiphcity. The 
two-particle methods measure V2,pp along the participant plane axis, except when 
the detector is sensitive to spectator neutrons in the reaction plane. 

Since nonflow effects are mainly due to few particle correlations, estimates of 
flow coefficients based on multi-particle correlations have the distinct advantage 
of reducing nonflow correlations. Higher order cumulants are multi-particle corre- 
lations where the contributions of lower order multiplets have been subtracted. In 
the cumulant method it can be shown ifTSl 1211 l22l |23]| that, for example, the four- 
particle correlation minus twice the square of the two-particle correlation eliminates 
two-particle nonflow effects : 

((m«,im«,2<,3<,4>> = {un.1Un.2K.3KA) - 2(m„,im* 2)^ = -v^{4} , (16) 

where the double brackets indicate the cumulant and m„ , is again the «''' harmonic 
unit vector of particle /. The fourth-root of this result is taken to get v„{4}. The 
statistical errors are larger than with the event plane method. A disadvantage is 
that can sometimes be negative, depending on the nature of flow fluctuations 
(see discussion of sensitivity to flow fluctuations in Sec. |2.4| i. An advantage is that 
the cumulant technique allows a consistent treatment of acceptance effects ifTSll . 
Normally, generating functions are used to calculate the cumulants |221 . but direct 
calculation is also possible [21 J . 

Multiparticle cumulants can also involve mixed harmonics 024II25L An important 
example here is the three-particle correlation 

(m«,im«,2mL,3> = ^«^2«, (17) 

which was successfully used at RHIC to suppress nonflow in the study of vi Il25ll26l 
and V4 1241 . 



2.3 q-distributions, Lee-Yang Zeros, Bessel and Fourier transforms 

The flow vector Q involves all the particles. In the absence of correlations, its length 
would grow as the square-root of the multiplicity M. Thus, to remove most of the 
multiplicity dependence, a reduced flow vector was defined II2TI as 

q„ = Q„/VM. (18) 

In the limit of M > 1 its magnitude is distributed HI] [15] |2T1 |22l as 

dN qn 2(72 / q„Vn\/M 

dqn \ 
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where Iq is a modified Bessel function. In this q-distribution method one looks at 
the length of the flow vector, not its angle |2|. The collective flow shifts the length 
distribution out by v,^ M and fluctuations broaden the distribution. Nonflow corre- 
lations reduce the effective multiplicity, thus also broadening the distribution. From 



just statistical effects, in Eq. 19 would equal 1 /2, but broadening the distribution 
increases CT„: 

(t2=^(1+M(t2J, (20) 

where 

al,^5n + 2ai. (21) 

The Lee-Yang Zeros method is an all-particle correlation designed to subtract 
nonflow effects to all orders 1231 |28l |29l . It is based on a 1952 proposal of Lee and 
Yang to detect a liquid-gas phase transition. Using the second-harmonic flow vector 
Q2, the projection on to an arbitrary laboratory angle d is 

M 

ef = ^w,cos[2(0,-0)], (22) 

!=1 

where the sum is taken over all the particles / with lab angles 0, and weights w,. 
Usually five equally spaced values of Q are used to average out detector acceptance 
effects. The essence of the method is to find a zero of a complex generating function, 
but in practice the first minimum of the modulus of the generating function along 
the imaginary axis is used. The sum generating function based on (which is a 
sum) is given by 

Gl{ir) = \{e"-Q')l (23) 

where r is a variable along the imaginary axis of the complex plane and the average 
is taken over all events. The square of the modulus is used to determine the first 
minimum. The position of the first minimum at the lab angle 9 is r^, and, for the 
case of unit weights, is related to the "integrated" flow by 

V^^m/rl (24) 

V2^lv^)jM, (25) 
where jqi = 2.405 is the first root of the Bessel function Jq and M is the multiphcity. 



In Eq. (25 1 the average is taken over the lab angles 0. As in all multi-particle meth- 
ods, this V2 is V2.RP 1301 [3TI . along the reaction plane axis. This can be thought of 
in this way: V2 is lower than (v2.pp) by the same amount as V2,rp is lower than V2,pp. 
A variant Lee- Yang Zeros sum generating function method has been devised l32l 
which produces an event plane Q-vector calculated with weights designed to elimi- 
nate autocorrelations and nonflow effects. 
The product generating function is 
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M 

Gf ('>) = I (Hll +'>^;Cos(2((/),- - 0))]) I . (26) 

.7=1 

Calculation of this generating function requires more computer time because the 
product over all particles has to be calculated for each value of r. Although the sum 
generating function works fine for V2, analyses for V4 (and vi |33|) relative to V2 
have to be based on the product generating function. This is because the product 
generating function is better at suppressing autocorrelation effects which are more 
important for mixed harmonics |28|. The Lee- Yang Zeros method only works for a 
sufficient signal-to-noise ratio. Since the signal is V2 and the noise is proportional 
to 1/a/M, the parameter % — V2VM determines the applicability of the method. It 
is found that the errors get large and the results scatter when % < 0.8, or the full 
event plane resolution is less than 0.6. When there is no flow the method will find 
a minimum from a fluctuation. For STAR at ^^s^ = 200 GeV AuH-Au 1231 the 
method fails for central collisions because V2 is small, and for peripheral collisions 
because the multiplicity is small. 

The method of Fourier and Bessel transforms 1 34 1 of the flow vector distributions 
is intimately related to the Lee- Young Zeros method; it clearly illustrates how the 
separation of nonflow effects happens. Let fi){Qn.x) denote the distribution in the x 
component of the flow vector (Eq. (|4])) for the case of zero flow, v„ = 0. Then, in the 
case of non-zero flow, and under the condition ^/M » 1, the corresponding distri- 
bution can be written as a superposition of /o distributions "shifted" in the direction 
of flow by an appropriate amount depending on the reaction plane angle IfTSI : 

dP f <if 

/(e„,,) = — — - / —MQn,x-VnMcos{n^)). (27) 
The Fourier transform of this distribution is: 

fik) = (e'^e-O = /^/ ^/G„,..e'^e"Yo(G«,..-v„Mcos(nf)) 

= / j dte'^'Mt) = h{kv„M)Uk). (28) 

Remarkably, the flow contribution is completely factored out, and the zeros of the 
Fourier transform are determined by the zeros of the Bessel function Jo{kv„M). One 
finds 

Vn=M/{kiM), (29) 

where ki is the first zero of the Fourier transform. The above result is the same as that 
obtained applying the Lee-Yang Zeros method using the sum generating function. 
In fact this relation to the Fourier transform was already pointed out in the original 
paper ES I. 

The two-dimensional Fourier transform of d^P/dQ„^dQ„ y 
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/■ dP 
= / dQnMkQn)— ^ Mkv„M), (30) 
J dQ„ 

is reduced to the Bessel transform of the distribution in the magnitude of the flow 
vector. Note that in this approach (valid in the Hmit of VM ^1) the flow contri- 
bution is decoupled from all other correlations, due to the collective nature of flow. 
Note also that in the same limit one expects the distribution of flow vectors to be 
Gaussian due to the Central Limit Theorem, thus explaining why fitting the distribu- 
tion to the form derived in Ref. ifTSl (such fits have been used in Ref. ||3][35]) is also 
not sensitive to nonflow correlations. Thus in this limit all three methods, the Bessel 
Transform, Lee-Yang Zeros, and fitting the ^-distribution, become very similar, if 
not equivalent. 



2.4 Methods comparison: sensitivity to nonflow and flow 
fluctuations 

The results obtained with different methods discussed in the previous section are 
affected by nonflow and flow fluctuations in different ways. Also, in some methods 
the results are closer to flow values in the participant plane and in others to the re- 
action plane. For example, correlating particles in the same rapidity range with a 
single harmonic would measure flow in the participant plane, but using a mixed har- 
monic method, with the first harmonic determined from spectator neutrons, would 
provide elliptic flow in the reaction plane. 

Nonflow 5„ is defined by Eq. These are correlations not associated with the 
reaction plane. Included in nonflow effects are jets, resonance decay, short-range 
correlations such as the Hanbury-Brown Twis (HBT) effect, and momentum con- 
servation 1361 l37l . There exists several methods to evaluate and suppress nonflow 
contributions, such as using rapidity gaps between correlated particles, using dif- 
ferent charge combinations for correlated particles (to assess the contribution of 
resonances), etc. 

Since nonflow correlations are mainly few-particle effects, 5„ roughly scales as 
the inverse of the multiplicit}[^ It leads to an almost constant contribution to the 
dependence of the expression M (uu*) on centrality while flow has a maximum for 
mid-central collisions, because in peripheral collisions the multiplicity is small and 
in central collision the anisotropic flow become small. Similarly, one can "subtract" 
nonflow contribution in flow measurement using the so called AA — pp method, as 
the nonflow contribution to the correlator (uQ*), is constant. For the scalar product 



This is true under assumption that the relative strengths of nonflow effects do not change with 
centrality. In reality one should probably expect some increase in nonflow effects in more central 
collisions due to larger relative contribution of hard parton collisions, and/or tighter azimuthal 
correlations from modification of the coirelation due to stronger radial flow 1 38 1. 
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method one can consider using 

{'^QlAA.correced = {^Q*)aA ~ i'^Ql pp , (31) 

where AA refers to nucleus-nucleus and pp refers to proton-proton. One can also 
assume a particular shape of nonflow correlations as a function of the difference 
in particle rapidity and/or azimuthal angles. Then one can estimate the flow value 
with the fit to the correlation function. We denote the results obtained from the fit to 
a 2-dimensional correlation function as V2{2d}. Note that although this notation is 
similar to one employed in Ref. |f39l . the meaning is different. 

Unfortunately, the above mentioned techniques do not allow real quantitative es- 
timates of the residual nonflow. Using multi-particle methods is more attractive in 
this sense, as they suppress nonflow effects by ^ \/M for each extra particle in 
the correlator. Estimates show that measuring elliptic flow at RHIC using 4-particle 
correlations almost completely removes nonflow effects. The largest remaining sys- 
tematic uncertainty is due to contributions to higher order cumulants from correla- 
tions when two particles, which are daughters of a resonance decay, are correlated 
with all other particles in the multiplet via flow of the resonance. Unfortunately, this 
effect can not be suppressed by using higher order cumulants. Global momentum 
conservation can affect measurements of directed flow, or elliptic flow measured 
with respect to the first harmonic event plane, if the detector acceptance is not sym- 
metric about mid-rapidity. It is a long-range effect and not reduced by a gap in 
pseudorapidity. The effect causes a discontinuity in vi at mid-rapidity, as seen in 
Fig. |5] If one can estimate the fraction of all produced particles which are detected, 
then a correction can be made for this effect ll40l . as is also shown in Fig. |5] Mo- 
mentum conservation is unimportant when the event plane is determined from an 
even Fourier harmonic, or for a detector having symmetric acceptance around mid- 
rapidity. 




Fig. 5 Directed flow as a function of rapidity for charged pions from minimum bias 158 A GeV 
Pb+Pb collisions l,40J . Shown are vi before (squares) and after (circles) correction for momentum 
conservation. 
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Flow fluctuations also affect different methods differently. The effect on cumu- 
lant results are the easiest to trace (though it does not mean it is easy to measure the 



fluctuations). The dependence follows directly from the definitions, Eqs. ( 14 16 1: 



v{2}=yV) = ((v)2 + (T2)V2, 

v{6} = ((l/4)((v«)-9(v^)(v2) + 12(v2)^))l/^ 



(32) 
(33) 
(34) 



etc. Note that although v'2{2} can be written in terms of only (v) and a^, ^{4} 
and higher cumulants in general require knowledge of higher order moments of the 
distribution in v. If needed, a model has to be used for the distribution of v to relate 
the contribution of flow fluctuations to different order cumulants. For example, for 
a Gaussian distribution |30| in v. 



V2{2} = ((v)2 + ct2)'/2;=.(v) + (t2/(2(v», 



V2{4} = ((v)^-2a2 



(v)-ct2/(2(v)), 



V2{6} = {{vr-3aUv)y^'^{v)-cy,'/{2{v)) 



(35) 
(36) 
(37) 



Note that the above relations are also valid for any other distribution in the limit 

C7v< (v). 

For v„ fluctuations according to a Bessel-Gaussian distribution: 



dN 



1 



the cumulants are ll30l: 



v^ + Vq 

' C7, 



72 



;BG(v„;vo,(J„) 



V2{2}2 = v2 + 2ct2, 
v'2{n} = v'o, n > 4. 



(38) 



(39) 
(40) 



Note, that the Gaussian model of elliptic flow fluctuations due to eccentricity fluc- 
tuations discussed in section |3.3.4| results in a Bessel-Gaussian distribution in V2 
with the parameter vq being V2.rp |30|, elliptic flow along the reaction plane axis. 



We also note that Eq. ( |T9] l is a Bessel-Gaussian for q, while Eq. (38 i is the Bessel- 
Gaussian for v„. Figures|6] and [7] are consistent with these equations based on the 
Bessel-Gaussian distribution. 

The event plane method dependence on flow fluctuations is more complicated 

1 /2 

and ranges from v'2{EP} — V2{2} = (v^) to V2{EP} = (v) depending on the re- 
action plane resolution |4T|. Defining parameter a via V2{EP} = (v" )'''" one finds 
that a « 2 for small values of resolution and approaches unity for large values. In 
Ref. BTl this observation was made based on Monte-Carlo simulations, but this 
dependence also can be obtained analytically in the case of small fluctuations, or 
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Fig. 6 Many particle cumulant results for 
charged hadron from 158A Pb+Pb as a func- 
tion of centrality with the most central at the 
left. Lines are polynomial fits 14). 
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Fig. 7 Many particle cumulant results for 
charged hadron V2 from y's^ = 130 GeV 
Au+Au as a function of centrality 1211 . 



numerically with direct integration over the v distribution in the numerator and de- 
nominator of Eq. ([8]l f42l- 

The dependence on fluctuations of the Lee- Yang Zeros and other similar meth- 
ods, such as fitting the ^-distribution or using Fourier-Bessel transforms, is also 
non-linear As was shown in Ref. |34| by Monte-Carlo simulations, the dependence 
is close to that of higher (n > 2) cumulants. In the small fluctuation limit it also can 
be obtained analytically. Note that for the Bessel-Gaussian distribution in v, these 
methods yield the same results as higher cumulants, namely (v) — vq of the Bessel- 
Gaussian. Thus, if the distribution of V2 is Bessel-Gaussian, all the multi-particle 
methods should give the same result: (v) = vq~ V2,rp 1301 . 

The above equations can be used for removing the "trivial" effect of flow fluctua- 
tions due to variation of elliptic flow within a wide centrality bin |21 1. As discussed 
below, real flow fluctuations 11211 l95l are very difficult to separate from nonflow 
effects. 

Methods comparisons are shown in Fig.|8]for charged hadron 200 GeV AuH-Au 
results. Plotted is the integrated V2 divided by the values for the standard event plane 
method v'2{EP}. This is an update of Fig. 29b from Ref. [^261 by the addition of val- 
ues for the Lee- Yang Zeros methods |23 1. The results fall into two bands: the upper 
band for the two-particle correlation results and the lower band for the multi-particle 
results. The event plane values are about 5% below most of the other two-particle 
results, and the pseudorapidity sub-event values for peripheral collisions drop be- 
low the other two-particle results. The pseudorapidity sub-event method apparently 
succeeds in removing more nonflow, especially for peripheral collisions. 

According to our current interpretation, the upper band represents averages along 
the participant plane with nonflow and fluctuation contributions, and the lower band 
averages along the reaction plane mostly free of nonflow and fluctuation contribu- 
tions. The separation of the two bands is a function of (7,o[: 



v2{2f-v2{'Xf = al, = 52 + 20^2- 



(41) 
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% Most Central 

Fig. 8 Charged hadron V2 divided by V2{EP}as a function of geometrical cross section for 
y/Sffff = 200 GeV Au+Au. The values are integrated over pj- for | TJ | < 1 .0, except for the 
LYZ product generating function which had | f] |< 1-3. Results are shown for the event plane 
method, random sub-events, pseudorapidity sub-events, scalar product, two-particle cumulants, 
four-particle cumulants, ^-distribution, and Lee-Yang Zeros sum generating and product generat- 
ing functions I23II26I . 



The event plane method is a special case with results being somewhere between 
V2{2} and V2{4} depending on the reaction plane resolution. For the higher order 
cumulants it can be seen in Figs. [6] and [7] that although the two-particle cumulant 
values are above the four-particle cumulants, all the still higher order cumulants 
agree with V2{4}. Figure |9] compares the four-particle cumulant, the Lee- Yang Ze- 
ros, and the event plane results as a function of pr- The cumulant values decrease 
with pt compared to V2{EP}as one would expect for a nonflow contribution from 
jets at high pr, but the Lee- Yang Zeros values do not for some unknown reason. 



3 Anisotropic flow: results and physics 
3.1 General 

It has been found that many physical processes contribute to the development of 
anisotropic flow during the evolution of the system created in the collision. Though 
different flow harmonics often reflect different physics, there are some common 
features, such as the mass dependence of the differential flow in the low transverse 
momentum region, or the role of coalescence, both of which we discuss in more 
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Fig. 9 A comparison of results for charged hadron four-particle cumulants and Lee- Yang Zeros 
with sum generating function, to the event plane method as a function of pj- for \ rj \ < 1.0 in 
10^0% most central Au-l-Au collisions at y^s^ = 200 GeV. The bottom panel shows the ratios of 
V2 divided by V2{EP} l23l . 



details later in this section. The entire collision is viewed as going through several 
stages: the formation of the initial state (the result of this stage is often referred as 
the "initial conditions") that takes time of the order of the spatial dimensions of the 
Lorentz contracted nuclei, thermalization period, (viscous) hydrodynamic expan- 
sion, "post-hydro" expansion that is most frequently simulated by a hadronic cas- 
cade, and finally chemical and kinetic freeze-out. These stages might not have well 
identified boundaries. The freeze-out stage is very likely continuous as a natural evo- 
lution of the cascade phase and not occurring suddenly. Formation of the anisotropic 
flow occurs during each and every one of these stages, but depending on the colli- 
sion energy, rapidity and transverse momentum, the importance of different stages 
changes. From the theoretical point of view our knowledge of different stages is also 
non uniform. 2d hydro is obviously much better studied than 3d. Though there has 
been significant progress in understanding the initial condition in the Color Glass 
Condensate (CGC) model, one would find that probably the least is known about 
the very first stages and in particular about the rapidity dependence of the "initial" 
conditions. Note that although several mechanisms could lead to very fast thermal- 
ization (at the level of a fraction of a fermi), the size of the nuclei after relativistic 
contraction ~ R/^Iihn / ^^/s^ is not negligible, and as noticed by Stock |43| could 
delay the start of elliptic flow development and consequently lower V2 values, except 
probably at the highest RHIC energies. 

Very little is known about flow fields at the end of the thermalization stage; this 
question is just starting to be explored. An example here would be the calculations 
of elliptic flow in the CGC approach, where the entire result might be thought of 
as due to the initial flow field. The results obtained in this model do not agree with 
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the data, but might be used as an initial condition for hydrodynamic calculations. 
Another interesting attempt l44ll (see also similar ideas in an earlier paper ||45l ) to 
obtain the initial flow field is via the so called "Landau matching condition" for 
the energy momentum tensor, assuming initial free streaming of partons with rapid 
thermalization at times of about one fermi. The possibility for such initial flow in 
particular allows one to start the hydrodynamic evolution at realistic times of about 
one fermi, and not loose the ability to reproduce large elliptic flow. 

In terms of initial longitudinal velocity fields, it might be important to take into 
account the initial velocity gradient along the impact parameter as illustrated in 



Fig. 10 Such a gradient directly contributes to the in-plane expansion rate (see 
Eq. 23 in Ref. |46|). This effect naturally also leads to directed flow (see the same 
Eq. 23), which is briefly addressed in |47 1. It will be very interesting to compare the 
calculations in such a model to very precise data from STAR P8ll on directed flow. 
The relation to other models B9l [50l predicting non-trivial dependence of directed 
flow on rapidity would be very interesting. Speculating on this subject one notices 
that viscous effects must also play an important role in such a scenario. 



Fig. 10 Initial longitudinal velocity profile in 
non-central nuclear collisions |46|. 




Fig. 11 Radial flow modulated by the elliptic 
component. The solid ellipse shows the overlap 
region and the arrows show the direction and 
magnitude of the expansion velocity. 



3.1.1 Interplay of anisotropic and radial flow 

In the collision scenario where the final particles are produced via freeze-out of the 
(locally) thermalized matter exhibiting collective flow, an interplay of radial expan- 
sion and anisotropic flow leads to a characteristic dependence of the differential 
flow v„{pt) on the mass of the particle ||5T][52|. At low transverse momenta heavier 
particles have lower v„{pt)', for large mass and low temperature the flow v„{pt) 
can become even negative (opposite to the integrated value). The mass "splitting" 
of the differential flow being dependent simultaneously on all three velocities (ra- 
dial, anisotropic and thermal) can be a sensitive test of the theoretical models, in 
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particular to the presence of the phase transition during the evolution of the system 
(see discussion below in section [J. 3. 

The physics of this phenomena is rather simple and depends on the relative mag- 
nitudes of the three velocities: average radial expansion velocity, anisotropic ve- 
locity (the amplitude of modulation in radial expansion velocity as function of the 
relative angle to the reaction plane), and thermal velocity, which depends on the 
temperature and the mass of the particle. In the case of elliptic flow, the anisotropic 
component is positive for the in-plane direction and negative for the out-of-plane 



direction, see Fig. 11 When thermal velocities are small, particles with low pr 
are produced mostly from regions of the source with low radial velocity, where the 
anisotropic flow velocity component is negative in order to partially compensate the 
radial flow. In the case when thermal velocities are large (light particles) compared 
to anisotropic velocity modulations, the effect becomes small. The description of 
the effect is most transparent in the blast wave model l35l ISTl l52l l53l . 

Blast wave fits to differential flow measurements of data have been first applied 
in Ref. ||54| . fitting proton v\{pt) in AuH-Au collisions at the AGS and reasonable 
values for anisotropic flow velocities were obtained. Elliptic flow was very success- 
fully fit in Ref. 1 35 1 where the model was further developed including a parameter 
responsible for the spatial geometry of the system. It is further discussed in sec- 
tion 3.3.6 The effect is stronger for heavier particles. It is noteworthy that elliptic 



flow of J/i/ |55 | and deuterons |56| become negative at low transverse momenta. 



3.1.2 Flow amplification by coalescence 

The distributions of particles produced via coalescence of some primordial particles 
should reflect the original distributions of the primordial particles. If the coales- 
cence process does not significantly affect the distributions of the original particles, 
then the standard coalescence formula can be applied (written here for the case of 
nucleon coalescence into light nuclei of atomic number A): 

E,d\,^fE^\'^ (42) 



where the coefficient Ba reflects the coalescence probability; in particular it includes 
the integration over the spatial distribution of nucleons and is inversely proportional 
to the correlation volume. For simpUcity in this equation it is assumed that protons 
and neutrons have similar momentum distributions. Neglecting for the moment pos- 
sible non-uniformity in the spatial distribution and concentrating on the momentum 
distribution, one finds that if nucleons are subject to anisotropic flow, the distribution 
in Eq. (|42]) of nuclei A becomes even more anisotropic. In this context coalescence 
leads to a simple scaling relation for flow: 



VnAPTA} ~AVn.p{pTAlA). 



(43) 
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The E877 Collaboration 1571 1581 at the AGS observed that the directed flow of 
deuterons indeed followed this coalescence rule. Scaling violations were observed at 
rapidity close to the beam rapidity, and attributed to the change in spatial distribution 
of nucleons in- and out-of-plane. Elliptic flow of deuterons was studied at RHIC by 
the PHENIX [|59j and STAR i56! Collaborations. Both collaborations found good 
consistency with predictions of a coalescence model. 

Anisotropic flow enhancement due to coalescence can be used also as a tool for 
the study of resonance regeneration, such as ot K*, during hadronic evolution of 
the system. The idea is that if regeneration is important, and a significant fraction of 
the resonances are produced via coalescence of other hadrons, then the elliptic flow 
of the resonances should be enhanced compared to that of direct production il60l . 

It appears that constituent quark coalescence |61 1 plays a very important role in 
particle production in the intermediate pj region; its relation to the formation of 



elliptic flow is discussed in section 3.3.7 



3.2 Directed flow 



3.2.1 Physics of directed flow 




Where the colliding nuclei start to overlap, dense matter is created which deflects 
the remaining incoming nuclear matter (see Figs. [T| and [T2]i. The deflection of the 
remnants of the incoming nucleus at positive rapidity is in the +x direction leading 
to (px) > 0, and the remnants of the nucleus at negative rapidity are deflected in 
the —X direction thus having a (p^) < 0. The magnitude of the deflection probes 
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Fig. 13 A sketch of a medium central symmetric heavy-ion collision progressing in time (a) and 
(c), and the rapidity distribution of (p^) and {x) in (d). In (b) the overlap region is magnified and 
the "spectators" are not shown. In these figures x is the coordinate along the impact parameter 
direction and z is the coordinate along the projectile direction J50J. 



the compressibility of the created dense matter. It probes the system at early time 
because the deflection happens during the passing time of the colliding heavy-ions. 
At AGS energies this is considered the dominant mechanism for generating directed 
flow. The observed directed flow at AGS and SPS energies is an almost linear func- 
tion of rapidity so that at these energies the slope dv'i {y)/dy at midrapidity was often 
used to quantify the strength of directed flow. 

At higher energies the linear dependence of directed flow is expected to break 
down; at midrapidity the directed flow is predicted to be very small and it is pos- 
sible that the slope at midrapidity has a sign opposite to that in the beam rapidity 
region. This so-called 'wiggle', whereby the directed flow changes sign three times 
outside the beam fragmentation region as illustrated in Fig. [TSjd), is very sensitive 
to the equation of state ISO] |62l |63l The wiggle can have different physical 
origins. Using a hydrodynamic approach it is observed in Refs. \ET^ "631 that this 
wiggle structure only appears under the assumption of a QGP equation of state, thus 
becoming a signature of the QGP phase transition. For another class of models, e.g. 



cascade models, the initial conditions are defined as illustrated in Fig. 13 Here the 
rapidity loss of the incoming nucleons at positive rapidity is larger at negative x than 
at positive x while for the incoming nucleons at negative rapidity the rapidity loss is 
larger at positive x than at negative x. This, in addition to a positive space momen- 
tum correlation, can also cause a wiggle structure in the directed flow [50 1. In this 
scenario the measured directed flow is sensitive to the magnitude of the rapidity loss 
and the strength of the space-momentum correlation, e.g. radial flow. More recently, 
taking into account this initial velocity gradient along the x-direction (see Fig.[TO]l, it 
is argued that one also has to consider the collective motion due to angular momen- 
tum conservation. It has been shown in Ref. f46l that this contributes to the in-plane 
expansion rate and as argued in Ref. [47J this should also lead to directed flow. 
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Fig. 14 Directed flow of charged hadrons for two beam energies and two colliding species as a 
function of pseudorapidity for 30-60% centrality. The event plane was determined from specta- 
tor neutrons in the shower maximumm detecors of the zero-degree calorimeters. The solid and 
dashed curves are odd-order polynomial fits to demonstrate the forward-backward symmetry of 
the measurements. The AMPT model calculations are plotted for the same conditions as the data 
(see legend, plotted only on one side of 7] = to reduce clutter). From Ref. f65i . 
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3.2.2 System size and energy dependence; extended longitudinal scaling 



Figure 14 shows the measured directed flow of charged hadrons at RHIC for two 
beam energies and two cofliding species. Clearly the magnitude of the directed flow 
is very small at midrapidity, as was expected, and the directed flow is not a linear 
function of rapidity as it was at lower beam energies. However there is also no sign 
of a wiggle structure in the observed charged particle directed flow. To rule out the 
models predicting such a wiggle, directed flow of identified particles has to be mea- 
sured because in the case of charged hadrons the wiggle could be masked due to the 
opposite sign of the wiggle for nucleons and pions as predicted in Ref. |50|. Fig- 
ure[T4]also clearly shows that at the same colliding energy and at the same fraction of 
cross section, the magnitude of directed flow is the same for Au+Au and Cu+Cu col- 
lisions over the whole measured pseudorapidity range. This is remarkable because 
the Au+Au system is three times more massive and indeed most model calculations, 
as shown for AMPT |661 in the figure, predict a stronger directed flow for the more 
massive system. In addition the figure shows that the directed flow decreases with 
collision energy at fixed rapidity. 



To investigate the collision energy dependence in Fig. 15 the directed flow is 
plotted as function of rj — ybeam- i e- in the frame of incoming nucleons. In this frame 
the directed flow for different species and collision energies is shown to collapse 
into a unified curve. The AMPT calculations clearly do not predict such a scaling; 
in AMPT the scaling only approximately holds for each colliding species separately. 
This remarkable scaling as a function of collision energy and for the two different 
colUding systems was to our knowledge not predicted in any model calculation. 
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Fig. 15 Directed flow of charged hadrons as a function of rj - 
and Cu+Cu collisions, each at two energies 1651. 
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However, the recent work in Refs. Il46ll47]| taking into account the initial flow fields 
might provide a plausible explanation. 



3.3 Elliptic flow 

3.3.1 In-plane elliptic flow 

Elliptic flow has attracted the most attention in recent years. Based on hydrodynamic 
calculations, in-plane elliptic flow was suggested fP] as a signature of collective ex- 
pansion in ultrarelativistic nuclear collisions. Note that the situation is different at 
lower energies, where the elliptical shape of the particle transverse momentum dis- 
tribution at midrapidity is elongated in the direction perpendicular to the reaction 
plane and interpreted as due to shadowing by spectator nucleons, a phenomenon 
called squeeze-out. Only at high energies, when the longitudinal size of the Lorentz 
contracted nuclei become negligible compared to the transverse size and, corre- 
spondingly, the passing time too small compared to the characteristic time for the 
development of elliptic flow, of the size of the nuclei radius, the shadowing goes 
away and elliptic flow fully develops in-plane. In-plane elliptic flow was observed 



at the AGS by the E877 Collaboration |3| (Fig.[T6]l and later, during the AGS low 
energy scan, the E895 Collaboration performed measurements of elliptic flow in 
the transition region from out-of -plane to in-plane. The first measurements of ellip- 



tic flow at RHIC 121] (Fig. 17 1 showed that V2 approached the predictions of ideal 



hydrodynamic models; this observation was taken as a signature of the rapid ther- 
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malization of the system. Figure [18] shows the excitation function of V2 including 
other lower beam energies and later measurements from SPS and RHIC. 




Fig. 16 First observation of in-plane elliptic flow of charged particles at 1 1.8 A GeV by the E877 
Collaboration 1 3 1 . The first harmonic event plane was determined by (mostly nucleon) directed 
flow detected by the Target Calorimeter. Mid-rapidity is 7] 1.7. 



" 0.1 




Fig. 17 First RHIC results [21] of elliptic flow plotted vs the charged hadron multiplicity divided 
by the maximum observed multiplicity. Blue boxes show estimates based on hydrodynamic calcu- 
lations. 
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Fig. 18 Elliptic flow for midrapidity 25% most central collisions as a function of beam energy 1671 . 



3.3.2 Low density and ideal hydro limits, v2/e plot 

As elliptic flow should be zero in the absence of rescattering, in the low den- 
sity limit (LDL) it is directly proportional to the particle density in the transverse 
plane ll68ll69l . As particle density increases, the anisotropy is expected also to in- 
crease, saturating at some level, assuming no modification to the constituent inter- 
actions. In addition, as the elliptic flow should be zero for a totally azimuthally 
symmetric system, for small anisotropies in the initial geometry, elliptic flow must 
be proportional to this spatial anisotropy. Usually, for this purpose one uses the ec- 
centricity defined by 

where the average is taken over the initial geometry with some weight. For the 
weight one can use the nuclear profile density (participant nucleons), entropy or 
energy densities, or something else (e.g. number of binary collisions). Note that 
as long as the eccentricity is small, elliptic flow should be directly proportional 
to eccentricity calculated with any weight, but it is important to keep weight the 
same when comparing results obtained in different calculations. For numerically 
large eccentricities the direct proportionality could break in principle, but as was 
shown in the very first hydrodynamic calculation by Ollitrault [ij | (though he used 
somewhat different measure of eccentricity) the proportionality holds well even for 
rather large values of e. It was pointed out by Sorge |,70J , who tried to study the 
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effect of the QGP phase transition based on RQMD calculations, that the centrality 
dependence of the scaled elliptic flow V2/£ (he used a different notation) would 
exhibit non-monotonic dependence in response to the softening of the equation of 
state. 

Based on all the above observations Voloshin and Poskanzer ll69l proposed to 
plot all the experimental data as v^js. vs particle density in the transverse plane, 
\ / S{dNch/ dy), where the initial overlap area S and eccentricity are taken from 
Glauber model calculations. The idea of the plot is to compare the results obtained 
at different collision energies, with different projectiles, and at different centralities. 
Non-smooth dependence would be indicative of new physics (for example, decon- 
finement) and saturation could signal an approach to ideal hydrodynamical evolu- 
tion. Figure[T9]shows NA49 results |4| together with results obtained at the AGS and 
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Fig. 19 Compilation of vj / e data ||4J vs particle density at midrapidity. Green lines indicate ideal 
hydrodynamic predictions for AGS, SPS and RHIC collisions energies 1711 . 



preliminary results from RHIC. The SPS and RHIC flow values were obtained with 
higher order cumulants and the eccentricities were taken from an optical Glauber 
model. As will be discussed in detail in section [3. 3. 4| this combination represents 
the best for comparison to theoretical calculations, as it is most free from both non- 
flow and fluctuations in the initial geometry of the overlap region. For this plot, 
elliptic flow values were integrated over the entire pj region; the data, if measured 
in a limited pj window were extrapolated to correct for this. Also, when rapidity 
density was not measured, the pseudorapidity density has been used with rescaling 
based on model calculations. For discussions of other systematic uncertainties see 
the original paper f4l . 

Figure 19 attracted a lot of attention as the data show a continuous rise reaching 
the ideal hydrodynamic expectations (shown by green lines) in the most central 
collisions at RHIC energies. The green arrow indicates the position of the color 
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percolation phase transition predicted by Satz |72|; unfortunately the systematic 
uncertainties in the experimental data (see discussion in Ref. ID) do not allow a 
definite statement whether the data exhibit any non-smooth behavior at this place. 
Future beam energy scans at RHIC should clarify this point. 

Strong elliptic flow observed in central Au+Au collisions at the highest RHIC 
energies, consistent with prediction of ideal hydrodynamics, lead to the picture of 
strongly coupled quark-gluon plasma, sQGP |73 |. Taking into account the signifi- 
cance of the V2/£ plot in helping to establish the sQGP picture, all components of 
this plot has been recently reevaluated and we now discuss in detail recent devel- 
opments in this area. Along with several indirect indications that in central AuH-Au 
collisions thermal equilibrium is not complete, it was found that very small viscous 
effects can lead to a significant reduction in the predicted elliptic flow compared 
to the ideal hydro case. There are significant effects even for viscosity equal to 
the conjectured lower limit of shear viscosity to entropy ratio (rj/s) iJ4\. Contin- 
uous freeze-out as implemented by a cascade afterburner also has a considerable 
effect IITSI . see Fig. 20 Lower values of elliptic flow would contradict experimental 



measurements if other effects, responsible for an increase of elliptic flow (compared 
to the "standard" hydrodynamic calculation) could not be identified. Several such 
effects have been reported: 

• First, it was demonstrated by Huovinen ll76ll that ideal hydro calculations, if tuned 
to describe spectra, yield larger elliptic flow than thought previously, which em- 
phasizes the need to describe the spectra and elliptic flow simultaneously. Fig- 
ure 21 shows calculations of elliptic flow in two scenarios; one (Partial Chemical 
Equilibrium) shown by the solid line is a fit to the spectra and significantly over- 
predicts elliptic flow. 

• It was shown, that in some models, e.g. CGC, the initial eccentricity can take 
significantly larger values than in the optical Glauber model that is usually used 
to set initial conditions in hydro calculations. The larger eccentricities inevitably 
lead to larger efliptic flow 1771 |78l |79l [80l . 

• Flow fluctuations, the nature of which is much better understood in recent years, 
lead to an increase of apparent flow relative to the participant plane |30|. 

• Finally, it was noticed that gradients in the initial velocity field (Fig. [T0| also 
increase the final values of elliptic flow B6l . 

Any of the above mentioned effects can be quite significant, each leading to 20- 
30% or even larger changes in the values of V2. The inclusion of these effects into 
one reliable model is still under way. 



3.3.3 Viscous effects 



An attempt of model independent analysis of V2 /e dependence on particle density 
based on a parameterization in terms of Knudsen number has been developed in 
Refs. ISTl |82I . By definition, l/K is, the mean number of collisions per particle, 
and ideal hydrodynamics corresponds to the limit \/K ^ oa. The authors use the 
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Fig. 20 Hydro+cascade calculations for v^je 
vs beam energy with CGC initial conditions 
and different freeze-out conditions. The points 
without a line through them were calculated 
with a cascade afterburner. Predictions for LHC 
are on the right side of the graph 1751 . 




Fig. 21 Ideal hydrodynamic calculations of 
pion i'2(pr) in two scenarios |76|. Solid line 
(Partial Chemical Equilibrium) indicates the 
results using parameters best fit to spectra. 
Dashed line assuming chemical equilibrium. 



expression V2/e = (i^'2 /e) hydro (1 +^^/-Ko) , where the parameter « 0.7 is inde 



pendently estimated from a fit of model calculations to the data (see Fig. 22 1. The 
authors conclude that at RHIC we might be up to 30% below the ideal "hydro limit" 
even for the most central collisions. Their estimate of the viscosity yields values 
of 7] /s = 0.11 —0.19 depending on the CGC or Glauber initial conditions. (These 
initial conditions are discussed more in Sec. |3.3.4| ) Similar fits to STAR data 
lead to similar conclusions. 



0.3 
0.25 

0.2 
0,15 

0.1 
0,05 




Au-Au 



CU-CU ' E3- 

fit a =4.3 mb 



...k 



,*-+■++■}■■ 



0.2 0.4 0.6 O.S 

(l/S)(llN..yy)[mb'- 



1.2 1.4 



Fig. 22 Fit to V2/e vs particle density in terms of Knudsen number for Cu-fCu and Au+Au L82.I 
The data are from PHOBOS using the eta subevent method. 



The magnitude of the viscous effects could be judged already from the early cal- 
culations ll84l where the hydrodynamical evolution at some intermediate stage was 
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joined to the transport model RQMDv2.4 to simulate late (viscous) evolution of 
the system and differential freeze-out. Using the transport code to describe the late 
hadronic stages of the system evolution and freeze-out, allowed a satisfactory de- 
scription of elliptic flow for SPS and RHIC data using the same EoS (QGP, including 
first order phase transition with a latent heat of 0.8 GeV/fm^), see Fig. 23 At SPS, V2 
values have been found about a factor of 2 lower compared to ideal hydrodynamic 
predictions. The importance of continuous freeze-out and late hadronic viscosity 



have been also demonstrated in similar calculations by Hirano, Fig. 20 Somewhat 



smaller in magnitude, but again a similar effect, was obtained in a hydroH-uRQMD 
hybrid calculation in Ref. ||85l . 

Taken together, hybrid model results show that in this approach one achieves 
very reasonable, and at the moment probably the best description of the data. The 
collision energy and centrality dependence in the hybrid models are mostly due to 
changes in the relative time the system spends in the sQGP state compared to the 
hadronic gas and the (continuous) freeze-out. Note that these models employ ideal 
hydrodynamics, once again arguing for very small viscosity of the early QGP stage. 
Viscose effects in this approach are totally due to the hadronic cascade phase. 
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Fig. 23 Hydro+RQMD results 1 86 1 for elliptic flow vs rapidity density. Experimental data are the 
two points at the top. LH8 denotes the results obtained for an EoS with latent heat 0.8 GeV/fm-'. 



Later, viscosity was attempted to be introduced directly into hydrodynamic cal- 
culations |87|. Recently there have been performed calculations 1^81 l90l 1911 l92l of 
the hydrodynamical expansion with viscous terms explicitly included into the equa- 
tions. Now everyone agrees ||93] on the significance of the viscous effects even for 
the conjectured minimal value of shear viscosity to entropy ratio (rj /s = 1/ (47r)). 
The results presented in Figs. 24 and 25 show that even minimal viscosity leads to 
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~25-30% reduction in flow values in Au+Au collisions and probably more than 
50% in Cu+Cu. 



♦ STAR EP 

ideal 

— r\/s=Om 
— ri''s=0.08 ; 

■ STAR cum/LYZ 



Fig. 24 Viscous hydro calculations [88' '891 compared to minimum bias STAR data |23| . The 
lower STAR points have an estimated correction for nonflow and fluctuation effects based on the 
four-particle cumulant and Lee- Yang Zeros methods shown in Fig. [9] 



Note that viscosity coefficients calculated in pQCD are usually much larger than 
would be allowed by the data. In this sense, noteworthy are the recent calcula- 
tions ||94ll which emphasize the importance of taking into account 2 3 processes. 
With these effects included, the viscosity coefficient appears to be about an order of 
magnitude smaller compared to previous calculations and falls into the "allowable" 
range of the data. 




Fig. 25 V2{pt) for 7t, K, and p particles from ideal and viscous (rj /s =1/ {47r)) hydrodynamics 
for Cu-l-Cu mid-central collisions |91| . 
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3.3.4 Initial eccentricity and V2 fluctuations 



The role of fluctuations in the initial system geometry f961 defined by nuclear par- 
ticipants 1.97] (interacting nucleons or quarks) has been greatly clarified in recent 
years ll30l ISTl l34l |4T1 l98l l99l . The following picture emerges: at fixed impact pa- 
rameter, the geometry of the participant zone (see Fig.|2| fluctuates, both, in terms 
of the value of the eccentricity as well as the orientation of the major axes. Then 
the anisotropy develops along the plane spanned by the minor axis of the participant 
zone and the beam direction, the so called participant plane. As the true reaction 
plane is not known and the event plane is estimated from the particle azimuthal 
distribution, the apparent (participant plane) flow appears to be always bigger (and 
always in-plane with V2.pp > 0) compared to the flow as projected onto the reaction 



plane V2,rp in Fig- 26 The importance of using the proper values of the initial ec- 



centricity is illustrated in Fig. 27 The V2 values from CuH-Cu and AuH-Au collisions 



are vastly different when scaled by Cstd (Eq. ( 44 1) but agree nicely when scaled by 



Epart (Eq.(45i)|[T00l. 
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Fig. 27 Left: Measured V2 of charged hadrons scaled by £std from Cu+Cu (upper curve) and 
Au+Au (lower curve) vs the number of participants. Right: Reseated with IIOOI . The data 
were obtained with the eta subevent method. 
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It was noticed in Ref. Il30l that in collisions of heavy nuclei the fluctuations in 
the eccentricity can be well described by a two-dimensional Gaussian with the same 
width in both directions: 



(^-%)= . h (45) 




where = (-*^) ~ i^)^- = (3'^) ~ (3')^' ^ty = (■^^3') ~ i^) iv)- What is not 
trivial is that for such Gaussian fluctuations the higher order cumulant flow values 
(v{n},n > 4) are not only insensitive to nonflow but also to eccentricity fluctuations. 
All of the higher cumulants are exactly equal to V2.rp, the flow projected onto the 
reaction plane. At the same time, the participant plane flow become unmeasurable 
because flow fluctuations can not be separated from nonflow contributions by means 
of correlation measurements. See Ref. P2]| . 

The fact that higher order cumulants of eccentricity are very close to the "stan- 
dard" values Estd was initially observed in numerical calculations |34|, and as an ap- 
proximate result in Ref. ||3TI . Numerically, deviations were observed only for lighter 
(CuH-Cu) systems, which were traced to a break in the Gaussian approximation |4l]. 
Experimentally, this question can be addressed by comparing the values of higher 
order cumulant flow with results obtained with the first harmonic reaction plane 
determined by spectator neutrons (at RHIC with the help of Zero Degree Calorime- 
ters). Note also that using the first harmonic reaction plane one can also address the 
validity of the Gaussian approximation by direct fluctuation measurements IIIOII . 

An important conclusion from that study of eccentricity fluctuations was that in 
most cases the measurements of higher order cumulant flow values provide elliptic 
flow relative to the reaction plane and not the participant plane. Similarly, the same 
value is given by several other methods, such as Lee- Yang Zeros, Bessel Transform, 
and the fit to q-distribution. This greatly simplifies the comparison of hydrodynamic 
calculations to data, as it says that in such calculations one should not worry about 
how to take into account fluctuations in the initial eccentricity (which is a non- 
trivial task) but just compare to the "right" measurement, V2,rp, e.g. V2{4}. This 
understanding also explains some earlier calculations with the uRQMD model II 1021 
[1031. There, it was shown that using higher cumulants and/or the LYZ method one 
indeed can measure elliptic flow very well, but it was not at all clear why there was 
no trace of the effects of flow fluctuations, which were expected in this model. 

Unfortunately this progress in understanding the nature of fluctuations does not 
help in resolving the problem of measuring flow fluctuations and nonflow. Strictly 
speaking, to separate those effects one is required to make assumptions. Most often, 
the azimuthal correlations between particles with large rapidity separation are used 
to suppress nonflow contributions. The problem with this method is that there are 
no reliable estimates of how well they suppresses nonflow and also how much the 
flow fluctuations (in this case correlations) change after imposing such a cut. 

At the QM'08 conference the PHOBOS lfT04l and STAR [105] coUaborations 
have presented their revised (compared to QM'06) results on flow fluctuations. 
These results are in good agreement, see Figs. 28 29 In Ref. 1 1051 a conservative 
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approach is taken and only upper limits on fluctuations are reported. The PHOBOS 
Collaboration uses estimates of nonflow effects from correlations with large rapid- 
ity separation and reports a more restrictive range for fluctuations. Both agree that 
the current measurements exhaust the (nucleon) eccentricity values obtained in the 
MC Glauber model and in this sense somewhat favor models which predict smaller 
relative fluctuations, such as the CGC model or a MC Glauber model taking into 
account constituent quark substructure. 



r,, T <'.,,./<e„^MC Glauber 

j\ nucleon participants 
2 B _ —*— canstiiuentqisarti pariicipants 
color glass (fKLN) 




impact parameter (b) (fm) 



Fig. 28 The shaded area shows the elliptic flow 
fluctuations as estimated by STAR 1 105|. The 
curves ai'e MC Glauber calculations for nucle- 
ons, quarks, and the CGC. 
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Fig. 29 The shaded area shows the eUiptic flow 
fluctuations as estimated by PHOBOS U04J . 
The Unes are estiinates from MC Glauber cal- 
culations. 



3.3.5 (Pseudo)rapidity dependence 

The measured pseudorapidity dependence of elliptic flow is shown in Figs. [30f[33| 
In general V2 exhibits a dNch/di] dependence, with a maximum at mid-rapidity and 
falls off in the fragmentation regions. Such behavior would find a natural explana- 



tion in the LDL. Collision energy and nuclear size dependence of V2{'r\), Fig. 30 



would also roughly agree with scaling proportional to the charged particle density. 



For 10^0% centrality Fig. 31 compares the PHOBOS and STAR results. It appears 
that using the Lee-Yang Zeros method has about the same effect as having a pseu- 
dorapidity gap. 

Similar to the limiting fragmentation picture of charged particle density, ellip- 



tic flow exhibits the same behavior. Figure 32 shows how the PHOBOS v'2(rj) data 
plotted against 77 — jbeam^ where jbeam is the projectile rapidity, collapses on the uni- 
versal curve in the fragmentation region. The original hydrodynamic models pre- 



dictions for V2{f]) are far from the data (see the solid curve in Fig. 33 1. Note that 
3d hydrodynamic calculations are computationally difficult and the dependence of 
the initial conditions on rapidity is far from being known. In all calculations so far a 
boost-invariant initial geometry is assumed. Adding a cascade hadronic afterburner 
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Fig. 30 Elliptic flow vs pseudorapidity at dif- 
ferent centralities 11061 . Insert shows ratio of 
peripheral to central. The data were obtained 
with the eta subevent method. 




Fig. 31 Pseudorapidity dependence of elliptic 
flow for y^s^ = 200 GeV Au-fAu from the 
event plane method (STAR), the LYZ method 
(STAR), and using an T] gap (PHOBOS). From 
Ref. IBl. 



greatly improves agreement with data (see Fig. 33 1 similar to the integral flow dis- 
cussed above (Figs . [20| and [23]l . 




Fig. 32 i'2 as a function of TJ — j'beam at four 
collisions energies |T00|. The data were ob- 
tained with the eta subevent methiod. 
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Fig. 33 Pseudorapidity dependence of ellip- 
tic flow from PHOBOS compared to a hy- 
dro-Fcascade model Ii75ll . The lines are for 
sudden freeze-out at different temperatures. 
The data were obtained with the eta subevent 
method. 
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3.3.6 Low pt region: mass splitting 



For charged hadrons, shown in Fig. 24 the elliptic flow increases almost linearly 
as a function of pr reaching values of about 0.15 at large pr- At low transverse 
momenta, the dependence is well described by hydrodynamics. 
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Fig. 34 Comparison of V2(pr) dependence of pions and protons with hydrodynamic calculations 
for a hadron gas and also including a phase transition t52lll07l . 



As explained in section 3.1.1 in a (locally) thermalized system, like in hydro- 
dynamics, the interplay of radial expansion and anisotropic flow should lead to 
a specific dependence of the differential flow vxipr) on the mass of the parti- 
cle iBTl l52l [T08J . Figure |34] shows V2 as function of transverse momentum for two 
particle species. As expected, at low pr the elliptic flow clearly depends on the 
mass of the particle with V2 at a fixed pj decreasing with increasing mass. The hy- 
drodynamic model calculations of V2{pt) for pions and (anti-)protons in Fig. 34 are 



performed for two equations of state: the full curves are for an EoS which incorpo- 
rates the effect of a phase transition from a QGP to a hadron gas, the dashed curves 
are for a hadronic EoS without phase transition. The hydro calculations clearly pre- 
dict the observed behavior rather well with a better description of the measurements 
provided by the EoS incorporating a phase transition. For the pions the effect of a 
phase transition is less pronounced compared to the protons. The lighter particles 
are more affected by the temperature, thus less sensitive to the collective flow veloc- 
ity and vice versa for the heavier particles. One should not, however, conclude from 
the good fit of the ideal hydro calculations to the data in Fig. 34 that the EoS which 
includes the phase transition is the only allowed EoS. To draw conclusions about 
the EoS one first has to better understand the initial conditions |[77l[78ll79ll8()l . has 
to have a more realistic description of the phase transition [109], and quantify the 
effects of viscous corrections ll87l [88l [89l |9U1 1911 l92l . 

Figure [35] shows that the elliptic flow of the different mass particles at low pr 
can be characterized rather well by a common set of four freeze-out parameters: the 
temperature, the mean radial flow velocity, the azimuthal dependence of the radial 
flow velocity and the source deformation ll35l . In hydrodynamics, these parameters 



Collective phenomena in non-central nuclear collisions October 30, 2008 Draft 37 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Pt [GeV/c] 

Fig. 35 Comparison of minimum bias vjipr) dependence on particle mass with blast-wave model 
fits taken from Ref. 1 110|. The data are a compilation of PHENIX and STAR V2{pt) measurements 
and the curves are from the fit to pion and proton V2 {pr ) performed in Ref. 1 35 1 . 



are not independent since they are related by the initial conditions and the equation 
of state. This also explains why V2{pt} differential flow for heavier particles is 
sensitive to the EoS. 

The four parameter blast-wave model fit, to the extent that it provides an accu- 
rate description of the system at freeze-out, has an advantage over hydrodynamic 
and cascade models in that it allows extraction of these freeze-out parameters with- 
out knowledge of the systems history before freeze-out. The blast-wave parameters 
obtained by fitting the pion and (anti-)proton V2{pt) in Ref- i351 also described the 
later obtained measurements of V2{pt) for the heavier A and Cascade particles, as 
can be seen in Fig. |35] In addition, the obtained S2 parameter which describes the 
source deformation was later confirmed also by azimuthally sensitive femtoscopy 
measurements [111 J. 

In ideal hydrodynamics the mass ordering in V2 persists up to large pr, although 
less pronounced because the V2 of the different particles start to approach each other. 
It is seen that at higher pr the measurements start to deviate significantly from hy- 
drodynamics for all particle species, and that the observed V2 of the heavier baryons 
is larger than that of the lighter mesons. This mass dependence is the reverse of the 
behavior observed at low pr- This is not expected for hadrons in hydrodynamics 
and is also not expected if the V2 is caused by parton energy loss (in the latter case 
there would be, to first order, no dependence on particle type). An elegant explana- 
tion of the unexpected particle type dependence and magnitude of V2 at large pr is 
provided by the coalescence picture 161111121 and is discussed in Sec. 3.3.7 
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3.3.7 Constituent quark number scaling 

Elliptic flow of identified particles measured in Au+Au collisions at RHIC exhibits a 
remarkable scaling with the number of constituent quarks - an apparent dependence 
of hadron elliptic flow at intermediate transverse momenta, pj ^2~A GeV/c, on 
the number of constituent quarks in the hadron. This observation is of particular 
interest and importance, as it indicates that the system is in a deconfined stage. In a 
more general sense, it appears that high energy nuclear collisions provide a window 
of opportunity to prove that hadron production indeed happens via the constituent 
quark phase. It has been noticed in Ref. |61| that if hadrons are formed via coa- 
lescence of the constituent quarks then there should be a region in the transverse 
momentum space where particle yield would be proportional to the quark density 
to the power equal to the number of constituent quarks in the produced hadron, 
2 for mesons and 3 for baryons. Besides other important consequences, such as 
enhanced relative production of baryons in this transverse momentum region, this 
picture leads to the constituent quark scaling of elliptic flow, V2{pt) ~ « vjipr/n), 
where n is the number of constituent quarks in the hadron iSTllI 121 . Figure[36]shows 
this scaling holds with good accuracy. Note that while the scaling itself is limited 
to a specific region in transverse momentum, the coalescence mechanism can be 
valid at all smaller momenta. The reason for the scaling violation at lower momenta 
might be that the equations used to describe coalescence break down. For identified 
particles the scaling for mesons and baryons is based on the equations 

which are valid only if the probability of coalescence is relatively low. At low trans- 
verse momentum, where most of the quarks hadronize via coalescence these equa- 
tions break unitarity. Note that according to these equations, the hadron yield scales 
with power 2 or 3 of the quark density. For a discussion of the particle spectra cal- 
culations and further development of the coalescence picture we refer to the review 
in Ref. fun. 

In this picture the quantity V2{pT/n) is interpreted as the elliptic flow of con- 
stituent quarks. This means deconfinement - as the constituent quarks must be in 
a deconfined phase in order to be freely "reshuffled" into final hadrons. This could 
be the first, and very strong argument for an observation of deconfined matter at 
RHIC. The meson is an interesting test of constituent quark scaling. It has the 



mass of a proton but contains only two quarks. Figure 37 shows that its V2 values 



are consistent with the two-quark curve thereby supporting the quark picture. 



Figure 38 shows that plotting V2 vs transverse kinetic energy KEj ~ nij — m. 



where mj = \J Pj + m^, results in the formation of distinct branches for mesons 
and baryons. Scaling with the number of constituent quarks n^ then coalesces the 
two branches into one curve. At the moment there is no agreement on the reason 
for such a universal scaling except that re-plotting the data vs transverse kinetic en- 
ergy to some extent compensates for the effect of radial flow (see the discussion 
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Fig. 36 Test of the constituent quark num- 
ber scaling of elliptic flow for minimum bias 
Au-l-Au collisions at ^s^ = 200 GeV. The 
dashed lines are polynomial fits Ql4J. 
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Fig. 37 Minimum bias V2 as a function of pr 
for the (p meson, K'^, and A + A fTTSl . The 
dashed and dotted lines represent what is ex- 
pected from coalescence for two- and three- 
quark particles. 
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Fig. 38 V2 for y's^ = 200 GeV Au-l-Au as a function of and KEj-, and also V2/n^ vs KEj- /n^ 
in the last panel 11161 . 



in section 3.1.1 1. But we mention that transverse kinetic energy scaling is an ap- 
proximation to a more general scaling found in the Buda-Lund model Ml 171 . Also, 
a quark recombination model based on q-qbar resonance interactions may account 
for this phenomenon ! 1 1 811 . 

A detailed study of necessary conditions for the quark-number scaling to be valid 
has been performed by Pratt and Pal 1 1 19 1. Two limiting possibilities are illustrated 
in Fig. [39] where the particles (quarks) of a given velocity are represented by arrows. 
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In case (a) the effective volumes of the right-moving and upward-moving particles 
are the same but the right-moving particles have a larger phase space density. In 
case (b) the densities are the same but the volumes differ. Both cases correspond to 
the same quark V2, but only case (a) would lead to the quark-number scaling. Based 
on this analysis one might conclude that constituent quark scaling contradicts local 
thermalization and freeze-out at a constant phase-space density. Even if true, this 
does not diminish the validity of the conclusion on deconfinement, but points to a 
very interesting possibility that the system created in the heavy ion collision can be 
in a deconfined but not completely thermalized state. It also does not exclude the 
possibility that thermalization happens only for lower transverse momenta. 




Fig. 39 Two possible source configurations, 
corresponding to the same quai'k elliptic flow, 
but only case (a) leads to quark-number scal- 
ing Gil)- 
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Fig. 40 Centrality dependence of elliptic flow 
in Au+Au collisions at y's^ = 130 GeV |26| 
compared to results of the AMPT model with 
string melting and different cross-sections 1, 66.1 . 



The success of quark coalescence as a hadronization mechanism in describing 
particle production at intermediate transverse momenta hints on the importance of 
the constituent quark dynamics in the evolution of the system in general. Recall that 
the transport models 1661 1 12011 in their standard configuration fail to describe the 
strong increase of elliptic flow with energy. They have to significantly increase the 
parton transport cross section or the matter density in order to reach the experimen- 
tal values. On the other hand if one looks carefully into what particular parameters 
are required in order to describe the data, an interesting picture emerges: the density 
and partonic cross sections are just what one would expect for a system of con- 
stituent quarks. For example. Fig. 40 shows a comparison of the experimental data 
to the AMPT |66 | model calculations in the so-called melted string scenario. The 
main assumption of this scenario is that the total number of partons in the system 
equals the number of constituent quarks in the produced hadrons, exactly what one 
would use for a model based on the picture of a system of constituent quarks. Note 
that the model describes the data best with a cross section of about 5 mb, which 
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is what one expects for constituent quarks. Similar conclusion can be drawn from 



Fig. 41 which shows the results of calculations in the MPC B120I model, using as 
a reference the standard HIJING gluon density (dNg/dy = 1200 for central Au+Au 
collisions) and a cross section of about 2 mb. Again, rescaling the parton density to 
that corresponding to the total number of quarks in the final hadrons and increasing 
the cross section to 5-6 mb one would achieve a reasonable description of the data. 




Fig. 41 Elliptic flow from MPC calculations compared to data for =130 GeV Au+Au fl20l|. 
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Fig. 42 Calculations 1 1 2 1 1 compared to experimental data. The LDL model results are shown by 
the dash-dot line. Other lines are hydrodynamic calculations with different equations of state. 
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Figure 42 shows the comparison of the experimental data on V2{pt) to calcula- 
tions in the low density limit II121L According to the equations of the LDL, V2{pt) 
saturates at transverse momenta about a few times the particle mass. The calcula- 
tions, shown by a dashed curve, exhibits this very behavior as first exhausting flow 
of pions and then at higher pr that of protons. Indeed the results do not resem- 
ble the data. But it is clear that similar calculations performed with masses of the 
constituent quarks (a few times that of the pion) taken together with coalescence 
enhancement at intermediate transverse momenta could well be consistent with the 
data. (See also the upper panel in Fig. 36 and the right panel of Fig. 38 which can 
be considered as V2{pt) of constituent quarks.) This again strongly advocates for a 
major role of constituent quark dynamics in the system evolution. 

Note that, in principle, hadronization always, even in small systems and at lower 
energies, occurs via constituent quark coalescence. The importance of the observa- 
tion of quark number scaling in elliptic flow, is that it constitutes deconfinement - 
"free" kinetics of quarks in the system. In this sense this scaling is a signature of 
deconfinement and should be violated under conditions when deconfinement does 
not happen. 



3.3.8 High pt region 

Elliptic anisotropy at high transverse momenta is an interesting observable as it 
is believed that it reflects the path length dependence of high pr parton energy 
loss 11221 . At sufficiently high transverse momentum in AuH-Au collisions, hadron 
yields are thought to contain a significant fraction originating from the fragmen- 
tation of high energy partons, resulting from initial hard scatterings. Calculations 
based on perturbative QCD (pQCD) predict that high energy partons traversing nu- 



clear matter lose energy through induced gluon radiation 11231 11241 . The energy 
loss (jet quenching) is expected to depend strongly on the color charge density of 
the created system and the traversed path length of the propagating parton. In non- 
central heavy-ion collisions, the geometrical overlap region has an almond shape 
in the transverse plane, with its short axis in the reaction plane. Depending on 
the azimuthal emission angle, partons traversing such a system, on average, ex- 
perience different path lengths and therefore different energy loss. This mechanism 
introduces an azimuthal anisotropy in particle production at high transverse mo- 
menta 02211 125! fT26J . 



The dashed Unes in Fig. 43 show the first quantitative theory predictions based on 
energy-loss calculations in a static medium [127 ] and are compared to STAR mea- 
surements 1 128 1 . However, as was discussed in the previous section, it was realized 
that in the pj range of 2 — 6 GeV/c hadron yields might not dominantly originate 
from the fragmentation of high energy partons but are instead mostly produced by 
quark coalescence. Therefore in order to compare to predictions of parton energy 
loss V2 has to be measured above pj = 6 GeV/c. Unfortunately nonflow contribu- 
tions which at low pr are modest are significant at high pr- In fact nonflow might 
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Fig. 43 The predicted magnitude of ^'2 at higher pj based on energy loss expectations for different 
values of the gluon density in a static medium |127|. The data were obtained with a modified 
event plane method where particles with 4 7] < 0.5 were excluded from the calculation of the event 
plane U28il . 



even start to dominate the measured two-particle azimuthal correlations at high pT 
which complicates making reliable V2 measurements. 
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Fig. 44 Azimuthal correlations in p+p, d+Au, and Au+Au collisions of different centralities 1261 
[1291. 



Assuming that there is no anisotropic flow in p+p and d+Au collisions one can 
estimate the contribution from nonflow by comparing the two- particle azimuthal 
correlation in these systems with the azimuthal correlation in Au+Au. Figure 44 
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shows the measured azimuthal correlation as a function of transverse momentum in 
p+p and d+Au colHsions as well as for three different centralities in Au+Au. The 
two-particle azimuthal correlation shown is defined in Eq. (31 1. It is seen that the 



magnitude of the azimuthal correlations in p+p, d+Au and very peripheral Au+Au 
is comparable over the complete measured transverse momentum range. This in- 
dicates that in peripheral Au+Au collisions nonflow correlations are the dominant 
contribution. On the other hand, for mid-central and central Au+Au collisions at 
low and intermediate pt the observed azimuthal correlation is much stronger than 
in p+p and d+Au as is expected in the presence of strong anisotropic flow. How- 
ever, at high Pt, even for these mid-central and central Au+Au collisions, the two- 
particle azimuthal correlation are again dominated by the correlations already seen 
in p+p and d+Au. The remaining difference in the azimuthal correlation at high pj 
can be used to estimate the elliptic flow, which is done in Refs. [26.,129J . 
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Fig. 45 Comparison of V2 (pr ) measured using 2- and 4- particle coirelations compai'ed to event 
plane results. The centrality was 20-60% for = 200 GeV Au+Au >1301 . 



Instead of using p+p and d+Au collisions as a reference for nonflow, the large 
data sample obtained at RHIC in AuH-Au collisions at y^s^ — 200 GeV now also 
allows for calculating the anisotropic flow coefficients up to 10 GeV/c using higher 



order cumulants. Figure 45 shows the V2{pt) obtained using the 4''' order cumulant, 
and the values are compared to the viipr) obtained using the event plane method 
and the 2-particle cumulant. The V2 values from the event plane method and 2- 
particle cumulant exhibit the same large values at high pr as seen in Fig. 44 The 
V2 obtained from the 4''' order cumulant method at high pT is at least a factor of 
two smaller, however still significant. These V2{^}{pt) values can be compared to 
model calculations to better constrain the mechanism responsible for parton energy 
loss. 

Another approach to reduce nonflow is to make sure that there is a large gap in 
rapidity between the particles used to determine the event plane and the particles 



Collective phenomena in non-central nuclear collisions October 30, 2008 Draft 



45 



I centralily 20%-60%~| 
> 




I centralily 20%-60% "| 



8 10 12 

P. (GeV/c) 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




- GuGu 900 GeV, STAR Preliminary 
: — ■ ^2(2} 

: • V2{FTPC} 




h 


















c* -I — n 






f 




■ 1 



8 10 12 

P. (GeV/c) 



Fig. 46 Comparison of V2 (pr ) measured from the 2-particle cumulant method and from the event 
plane method using the STAR FTPCs. The centrality was 20-60% for = 200 GeV 11311 . 

Au+Au is on the left and Cu+Cu on the right. 



used to calculate the elliptic flow. In STAR, in addition to the ZDCs, the Forward 
TPCs (FTPCs) which are located at 2.9 < Tj < 3.9 can be used for this 1131]. Fig- 
ure 46 left compares the V2{pt) measured at midrapidity (77 < |0.9|) from the 2- 



particle cumulant and from the event plane method using the STAR FTPCs. As was 
already observed, the nonflow (here the difference between V2{2} and V2{FTPC}) 
is significant above 2 GeV/c and becomes large above 6 GeV/c. Assuming that the 
nonflow is completely removed with Arj ^ 2 the difference between V2{4} from 
Fig. 45 and V2{FTPC} in Fig. 46 left shows the difference between V2 in the reac- 
tion plane and V2 in the participant plane with fluctuations. However, because the 
curves still turn up at very high pr, there is probably still some nonflow even with 
the event plane determined in the FTPCs. 

Clearly for mid-central collisions in Au+Au the nonflow contribution at high pj 
is already significant. In case of Cu+Cu, i.e. a smaller system, the elliptic flow for 
the same midcentral collisions is expected to be smaller while the relative nonflow 
contribution is expected to be larger Figure 46 right shows the V2{pt) measured 
in Cu+Cu with the 2"^^ order cumulant and with the event plane method using the 
STAR FTPCs fl311. It is seen that the V2 in Cu+Cu obtained from the two-paiticle 
cumulant method at intermediate pj (4-12 GeV/c) is even larger than in Au+Au. 
In contrast the V2{pt) in Cu+Cu obtained from the event plane method using the 
FTPCs is smaller than in Au+Au over almost the whole pj range. This confirms that 
in Cu+Cu indeed the V2{pt) is smaller than in Au+Au while the nonflow effects are 
larger. Even though statistical uncertainties are large for V2{FTPC} above 6 GeV/c 
the transition to again increasing values of V2 might indicate that a rapidity gap of 
4 T] ~ 2 is still not large enough to remove all nonflow in Cu+Cu. 
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3.3.9 Rare probes 

Elliptic flow of heavy flavor and direct photons, often called by the common name 
rare probes, attract significant attention as they provide important additional infor- 
mation about the dynamics and properties of the sQGP created in heavy ion col- 
lisions. Heavy flavor elliptic flow is studied either directly or via single electrons 
from semi-leptonic decays. 
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Fig. 47 Heavy flavor electron eUiptic flow L134il . The shaded area shows the model prediction 
based on a relativistic Langevin simulation of heavy quarks 11351 . 



Elliptic flow of hadrons containing charm or bottom quark(s) provide direct in- 
formation on the collectivity of heavy flavor quarks in the deconfined stage. So far 
there exist data from PHENIX on flow of single electrons and very limited data on 
flow of J/^. Simulations show that for pj > 2-3 GeV/c most electrons are from 
semi-leptonic decays of charm and bottom, and at pr > 5-8 GeV/c it could be dom- 
inated by bottom decays. At such transverse momentum, electrons to a large degree 
preserve the direction of heavy flavor mesons and thus can be studied to evaluate the 
elliptic flow of charm or bottom quarks. In a simple coalescence picture B132II the 
elliptic flow of charm hadrons is a combination of the elliptic flow of corresponding 
quarks taken at appropriate transverse momentum (roughly shared in proportion to 
the quark masses). It appearers that, for example, charm meson elliptic flow at trans- 
verse momentum of about 3 GeV/c may differ more than a factor of 2 depending on 
whether charm quarks participate in collective flow 11321 . 

PHENIX measurements III33II1341 . shown in Fig. 47 indicate that heavy quarks 
flow at the same level as light ones. From this measurement (combined with heavy 
flavor electron spectra) PHENIX derives an upper limit on viscosity rj/s < (1.3 — 
2)/(4;r), confirming that the system behaves as a near perfect liquid. Results on 
elliptic flow of J have just started to appear |134|, and in spite of large errors 
they immediately attract attention as they are largely consistent with out-of-plane 
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Fig. 48 V2(pr) of thermal photons 11371 in Au+Au collisions at ^/s^ = 200 GeV at = 0. 
Quark and hadronic matter contributions are shown separately. Elliptic flow of pions and p mesons 
obtained with the pair-wise method are shown for comparison. 
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Fig. 49 Photon V2{pt) for Au+Au collisions at RHIC fT79| obtained with the event plane de- 
tennined from the beam-beam counters at | T] | from 3.1 to 3.9. The dashed line shows jet- 
fragmentation and induced bremsstrahlung only, while the solid lines give jet-photon conversion, 
primary hard and thermal photons. The dotted line shows direct photons from decay of neutral 
mesons from jet fragmentation. The dot-dashed line also includes photons from pion decays. 



elliptic flow, thus possibly indicating strong effects of radial flow II136I . as discussed 
in Sec. imi 

Elliptic flow of direct photons has a very rich physics as there have been identi- 
fied several mechanisms which lead to significant elliptic flow. The most important 
are: (i) thermal photons, (ii) photons from high pj parton fragmentation (in vac- 
uum), (iii) bremsstrahlung and "jet conversion" photons from fast partons propagat- 
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ing inside the medium. Thermal photons originate both in QGP and hadronic phases. 
Calculations B137I show that even though the QGP phase is short due to higher tem- 
perature, it gives the largest contribution to photons at intermediate transverse mo- 
menta (see Fig. 48 1. High pr photons come from hard parton-parton scattering with 
subsequent usual fragmentation (leading to positive V2) or to bremsstrahlung, which 
is stronger out-of-plane leading to negative V2- The combination of all these effects 
results in a rather complicated picture (see Fig. 49 1. There is hope that with large 
statistics different mechanisms could be separated with different isolation cuts. Fi- 
nally, note that the mechanism, similar to that of "jet conversion" could also lead to 
an enrichment of jet fragments with some rare particles (including strangeness con- 
tent). It would lead to the corresponding negative contribution to the elliptic flow of 
that particle lfT38]| . 



3.4 Higher harmonics 



The higher harmonics give more detail about the shape of the azimuthal anisotropy. 
Values for v'4 and vg are shown in Fig 



50 



|24|. The values appear to scale as v„ oc y 



n/2 



as shown in Fig. 50 For identified particles in Fig. 51 V4 appears to have the same 



scaling as V2 with transverse kinetic energy KEj-, but vs KEr/n^ it seems to scale 
as the square of the number of constituent quarks II 1401 as one would guess because 
V4 °= V2 II24II . The proportionality constant can be seen in Fig. 52 For ideal hydro 
this ratio should be 1/2 [108] at high pj. Boltzmann calculations indicate that a 



finite Knudsen number 
non-viscous hydro. 



141J is needed to describe the deviation of the data from 
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Fig. 50 V2, V4, and \y, vs pr with respect to the second harmonic event plane, for minimum bias 
200 GeV Au+Au collisions. Also shown is V4 from three-particle cumulants. The dotted curves are 
1.2 vl and 1.2 v| Ell- 



in simple coalescence models II142I . the ratio V4 for hadrons is related to V4 
for quarks: 
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Fig. 51 V4 for 20-60% centrality y's^ = 200 GeV Au+Au as a function of pj- and KEj-, and 
also V4 / (ttcy ) ^ vs KEj- / in the last panel I140i . The event plane was determined with the RxNP 
detector at | r] | from 1.0 to 2.8. 
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Fig. 52 V4/V2 vs pr for 20-30% centrality for ^/sj^ = 200 GeV Au+Au. Also shown are ideal 
hydro and Boltzmann calculations with various Knudsen numbers I83II130I . 




h/vi]r;r-l/4+(l/2)[v4/v^]7\ (47) 
h/v^]3r"«l/3 + (l/3)[v4/vi]7\ (48) 

where here pr is the quark pr- The V4/V2 for mesons can also be related to V4/V2 
for baryons: 

hA'^]3r"-l/6+(2/3)[v4/vi]^;;". (49) 
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From the results for identified particle V4 it appears that the quark V4 / is approxi- 
mately 2 lfTT4ll . 



4 Conclusion and outlook 

The wealth of information obtained recently in the area of anisotropic flow is far 
from being fully explored. The theoretical understanding of the experimental data is 
rapidly developing and greatly contributes to the overall picture of heavy ion colli- 
sions dynamics as well as the properties of the new state of matter, sQGP. Unfortu- 
nately, we have not been able to discuss many other important measurements related 
to the event anisotropy, such as azimuthally sensitive identical and non-identical par- 
ticle femtoscopic measurements 11431 11441 (for a recent review, see Ref. II145II ) or 
two particle correlations relative to the reaction plane II 1291 1 146 1. 

The next years promise to bring new very important data which will shed more 
light on quite a few remaining questions. The most important developments we 
expect are in the following directions: (i) upgrade of RHIC detectors will give more 
precise information on identified particle anisotropic flow, (ii) beam energy scan at 
RHIC will allow looking very carefully in the region where critical phenomena are 
expected (e.g. color percolation, critical point), in particular to look for any non- 
monotonic behavior or scaling violation in the v^jE plot, (iii) Uh-U collisions at 
RHIC, (iv) and, of course, the first data at LHC is eagerly expected to answer what 
happens at higher beam energies. 

The main interest in the RHIC beam energy scan is the search for the QCD crit- 
ical point. The scan would cover the energy region from top AGS energies, over 
the CERN SPS energies, and higher. In terms of anisotropic flow the major observ- 
ables to watch would be a possible non-smooth behavior in the vijz dependence 
on particle density |69|, a possible disappearance of constituent quark number scal- 



ing, and the "collapse" of directed flow |147| (see Fig. 53 1. RHIC also has plans 
to extend its reach in terms of energy density using uranium beams. From the first 
estimates and ideas of using uranium beams we now have detailed simulations 11481 
of such collisions with methods developed for selection of the desired geometry of 
the collision. 

The predictions for the LHC are rather uncertain, though most researchers agree 
that elliptic flow will continue to increase |149|, partially due to the relatively 



smaller contribution of viscous effects. Simple extrapolations 1150111511 in Fig. 54 



of the collision energy dependence of V2 look rather reliable. Theoretical calcula- 
tions predict significantly smaller, thought still finite, viscous effect at LHC ener- 



gies. One can see it in Fig. 20 which shows V2/e evolution with particle density up 
to LHC energies in 3d-hydroH-cascade approach. Note that there exist some calcula- 
tions predicting a decrease of the elliptic flow 1 15^. Another important observation 
is an increase in mass dependence (splitting) of V2{pt) due to a strong increase of 
radial flow. 
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In summary, we have had very exciting years of anisotropic flow study, which 
greatly enriched our understanding of ultra-relativistic nuclear collisions and multi- 
particle production in general. We are looking forward to more new physics with 
RHIC and LHC. 




[AGeV] 

Fig. 53 "Collapse" of directed flow as dis- 
cussed in Ref. 11471 . The old measure of flow, 
the slope at mid-rapidity, is shown as a function 
of beam energy with an arrow line through the 
data. 
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Fig. 54 Extrapolation of V2 { J] } values to LHC 
energies 1 150| based on limiting fragmentation 
hypothesis and assuming linear dependence of 
elliptic flow on I ?7 — ybeam | ■ The x-axis is r] — 
>'beain shifted by the y of the LHC. The measure- 
ments are from the eta subevent method. 



5 Acknowledgments 

We thank all our colleagues for numerous discussions on many of the topics dis- 
cussed in this review. We thank our colleagues in STAR for permission to use some 
of our text from published STAR papers. We appreciate comments on the paper from 
Ante Bilandzic, Pasi Huovinen, Claude Pruneau, Ralf Rapp, and Gang Wang. This 
work was supported in part by FOM and NWO of the Netherlands, and the HENP 
Divisions of the Office of Science of the US Department of Energy under Contract 
Numbers DE-AC02-05CH1 1231 and DE-FG02-92ER40713. 



52 



S.A. Voloshin, A.M. Poskanzer, and R. Snellings 



References 



1. J. Y. Omtrault, Phys. Rev. D 46, 229 (1992). 

2. J. Barrette et al. [E877 Collaboration], Phys. Rev. Lett. 73, 2532 (1994) [arXiv:hep- 
ex/9405003]. 

3. J. Barrette et al. [E877 Collaboration], Phys. Rev. C 55, 1420 (1997) [Erratum-ibid. C 56, 
2336 (1997)] [arXiv:nucl-ex/96 100061. 

4. C. Alt et al. [NA49 Collaboration], Phys. Rev. C 68, 034903 (2003) [arXiv:nucl-ex/0303001]. 

5. K. H. Ackermann et al. [STAR Collaboration], Phys. Rev. Lett. 86, 402 (2001) [arXivmucl- 
ex/0009011]. 

6. K. G. Boreskov, A. B. Kaidalov and O. V. Kancheh, arXiv:0809.0625 [hep-ph]. 

7. B. Z. KopeUovich, A. H. Rezaeian and I. Schmidt, arXiv:0809.4327 [hep-ph]. 

8. D. Teaney and R. Venugopalan, Phys. Lett. B 539, 53 (2002) [arXiv:hep-ph/0203208]. 

9. A. Krasnitz, Y. Nara and R. Venugopalan, Phys. Lett. B 554, 21 (2003) [arXiv:hep- 
ph/0204361]. 

10. W. Reisdorf and H. G. Ritter, Ann. Rev. Nucl. Part. Sci. 47, 663 (1997). 

11. N. Herrmann, J. P. Wessels and T. Wienold, Ann. Rev. Nucl. Part. Sci. 49, 581 (1999). 

12. S. Voloshin and Y. Zhang, Z. Phys. C 70, 665 (1996) [arXiv:hep-ph/9407282]. 

13. J. Y. OUitrault, arXiv:nucl-ex/971 1003. 

14. J. Y. OUitrault, Nucl. Phys. A 638, 195 (1998) [arXiv:nucl-ex/9802005]. 

15. A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 1671 (1998) [arXiv:nucl-ex/9805001]. 

16. P. Danielewicz and G. Odjmiec, nucleusnucleus collisions" Phys. Lett. B 157, 146 (1985). 

17. J. Y. OUitrault, Phys. Rev. D 48, 1132 (1993) [arXiv:hep-ph/9303247]. 

18. N. Borghini, P. M. Dinh and J. Y. OUitrault, Phys. Rev. C 64, 054901 (2001) [arXiv:nucl- 
th/0105040]. 

19. I. Selyuzhenkov and S. Voloshin, Phys. Rev. C 77, 034904 (2008) [arXiv:0707.4672 [nucl- 
th]]. 

20. S. Wang et al, Phys. Rev. C 44, 1091 (1991). 

21. C. Adler et al. [STAR Collaboration], Phys. Rev. C 66, 034904 (2002) [arXivmucl- 

ex/0206001]. 

22. N. Borghini, P M. Dinh and J. Y. OUitrault, arXiv:nucl-ex/0110016. 

23. B. I. Abelev et al. [STAR CoUaboration], Phys. Rev. C 77, 054901 (2008) [arXiv:0801.3466 
[nucl-ex]]. 

24. J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92, 062301 (2004) [arXivmucl- 
ex/0310029]. 

25. N. Borghini, P M. Dinh and J. Y. OlUtrault, Phys. Rev. C 66, 014905 (2002) [arXivmucl- 

th/0204017]. 

26. J. Adams et al. [STAR CoUaboration], Phys. Rev. C 72, 014904 (2005) [arXivmucl- 

ex/04090331. 

27. J. Y. OUitrault, Nucl. Phys. A 590, 561C (1995). 

28. R. S. Bhalerao, N. Borghini and J. Y. OUitrault, Nucl. Phys. A 727, 373 (2003) [arXiv:nucl- 
th/0310016]. 

29. N. Borghini, R. S. Bhalerao and J. Y. OUitrault, J. Phys. G 30, S1213 (2004) [arXivmucl- 
th/0402053]. 

30. S. A. Voloshin, A. M. Poskanzer, A. Tang and G. Wang, Phys. Lett. B 659, 537 (2008) 
[arXiv:0708.0800 [nucl-th]]. 

31. R. S. Bhalerao and J. Y. OUitrault, Phys. Lett. B 641, 260 (2006) [arXiv:nucl-th/0607009]. 

32. A. Bilandzic, N. van der Kolk, J. Y. OUitrault and R. Snellings, arXiv:0801.3915 [nucl-ex]. 

33. N. Borghini and J. Y. OUitrault, Nucl. Phys. A 742, 130 (2004) [arXiv:nucl-th/0404087]. 

34. S. A. Voloshin, arXiv:nucl-th/0606022. 

35. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87, 182301 (2001) [arXivmucl- 
ex/0107003]. 

36. N. Borghini, R M. Dinh and J. Y. OlUtrault, Phys. Rev. C 62, 034902 (2000) [arXivmucl- 
th/0004026]. 



Collective phenomena in non-central nuclear collisions October 30, 2008 Draft 



53 



37. P. M. Dinh, N. Borghini and J. Y. OUitrault, Phys. Lett. B 477, 51 (2000) [arXiv:nucl- 
th/9912013]. 

38. S. A. Voloshin, Phys. Lett. B 632, 490 (2006) [arXiv:nucl-th/0312065]. 

39. T. A. Trainor and D. T. Kettler, Int. J. Mod. Phys. E 17 (2008) 1219 [arXiv:0704.1674 [hep- 
ph]]. 

40. N. Borghini, P. M. Dinh, J. Y. OUitrault, A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 
66, 014901 (2002) [arXiv:nucl-th/0202013]. 

41. B. Alver et at., Phys. Rev. C 77, 014906 (2008) [arXiv:071 1.3724 [nucl-ex]]. 

42. J.-Y. OUitrault, A.M. Poskanzer, S.A. Voloshin, in preparation, 2008. 

43. R. Stock, private communication, 2008. 

44. W. BroniowsM, M. Chojnacki, W. Florkowski and A. Kisiel, Phys. Rev. Lett. 101, 022301 
(2008) [arXiv:0801.4361 [nucl-th]]. 

45. M. Gyulassy, Yu. M. Sinyukov, I. Karpenko and A. V. Nazarenko, Braz. J. Phys. 37, 1031 

(2007). 

46. F. Becattini, F. Piccinmi and J. Rizzo, Phys. Rev. C 77, 024906 (2008) [arXiv:071 1.1253 
[nucl-th]]. 

47. S. M. Troshin and N. E. Tyurin, arXiv:0709.4090 [hep-ph]. 

48. G. Wang, J. Phys. G 34, S1093 (2007) [arXiv:nucl-ex/0701045]. 

49. L. P Csemai, J. I. Kapusta and L. D. McLerran, J. Phys. G 32, SI 15 (2006). 

50. R. J. M. SnelUngs, H. Sorge, S. A. Voloshin, F Q. Wang and N. Xu, Phys. Rev. Lett. 84, 2803 
(2000) |arXiv:nucl-ex/9908001]. 

51. S. A. Voloshin, Phys. Rev. C 55, 1630 (1997) |arXiv:nucl-th/961 1038]. 

52. P. Huovinen, P. F. Kolb, U. W. Heinz, R V. Ruuskanen and S. A. Voloshin, Phys. Lett. B 503, 
58 (2001) |arXiv:hep-ph/0101136]. 

53. F Retiere and M. A. Lisa, Phys. Rev. C 70, 044907 (2004) [arXiv:nucl-Ui/03 12024]. 

54. J. Barrette et al. [E877 Collaboration], Phys. Rev. C 56, 3254 (1997) [arXivmucl- 
ex/9707002]. 

55. A. Franz [PHENIX CoUaboration], arXiv:0805.1718 [nucl-ex]. 

56. H. d. Liu, J. Phys. G 34, S1087 (2007) [arXiv:nucl-ex/0701057]. 

57. S. Voloshin et al. [The E877 CoUaboration], Nucl. Phys. A 638, 455C (1998) [arXiv:nucl- 
ex/9802001]. 

58. J. Barrette et al. [E877 CoUaboration], Phys. Rev. C 59, 884 (1999) [arXiv:nucl-ex/9805006]. 

59. S. Afanasiev et al. [PHENK CoUaboration], Phys. Rev. Lett. 99, 052301 (2007) [arXivmucl- 

ex/0703024]. 

60. J. Adams et al. [STAR CoUaboration], Phys. Rev. C 71, 064902 (2005) [arXivmucl- 
ex/0412019]. 

61. S. A. Voloshin, Nucl. Phys. A 715, 379 (2003) [arXiv:nucl-ex/0210014]. 

62. J. Brachmann et al, Phys. Rev. C 61, 024909 (2000) [arXiv:nucl-th/9908010]. 

63. L. R Csemai and D. Rohrich, Phys. Lett. B 458, 454 (1999) [arXiv:nucl-th/9908034]. 

64. L. P Csernai et al.. Acta Phys. Hung. A 22, 181 (2005) farXiv:hep-ph/0405277]. 

65. B. 1. Abelev et al. [STAR Collaboration |, arXiv:0807.1518 [nucl-exj. 

66. Z. w. Lin and C. M. Ko, Phys. Rev. C 65, 034904 (2002) [arXiv:nucl-th/0108039]. 

67. A. Wetzler, private communication (2005). 

68. H. Heiselberg and A. M. Levy, Phys. Rev. C 59, 2716 (1999) [arXiv:nucl-th/98 12034]. 

69. S. A. Voloshin and A. M. Poskanzer, Phys. Lett. B 474, 27 (2000) [arXiv:nucl-th/9906075]. 

70. H. Sorge, Phys. Rev. Lett. 82, 2048 (1999) [arXiv:nucl-th/9812057]. 

71. P F Kolb, J. SoUfrank and U. W. Hemz, Phys. Rev. C 62, 054909 (2000) [arXiv:hep- 
ph/0006129]. 

72. H. Satz, Nucl. Phys. A 715, 3 (2003) [arXiv:hep-ph/0209181]. 

73. M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005) [arXiv:nucl-th/0405013]. 

74. P Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005) [arXiv:hep- 

th/0405231]. 

75. T. Hirano, Prog. Theor. Phys. Suppl. 168, 347 (2007) [arXiv:0704.1699 [nucl-th]]. 

76. R Huovinen, arXiv:07 10.4379 [nucl-th]. 



54 



S.A. Voloshin, A.M. Poskanzer, and R. Snellings 



77. T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey and Y. Nara, Phys. Lett. B 636, 299 (2006) 
[arXiv:nucl-th/0511046]. 

78. H. J. Drescher, A. Dumitm, A. Hayashigaki and Y. Nara, Phys. Rev. C 74, 044905 (2006) 
farXiv:nucl-th/06050121. 

79. H. J. Drescher and Y. Nara, Phys. Rev. C 76, 041903 (2007) [arXiv:0707.0249 [nucl-thl]. 

80. T. Lappi and R. Venugopalan, Phys. Rev. C 74, 054905 (2006) [arXiv:nucl-th/0609021 1. 

81. R. S. Bhalerao, J. P Blaizot, N. Borghini and J. Y. Ollitrauk, Phys. Lett. B 627, 49 (2005) 
[arXiv:nucl-th/0508009]. 

82. H. J. Drescher, A. Dumitru, C. Gombeaud and J. Y. OUitrault, Phys. Rev. C 76, 024905 
(2007) [arXiv:0704.3553 [nucl-th]]. 

83. A. Tang [for the STAR Collaboration], arXiv:0808.2144 [nucl-ex]. 

84. D. Teaney, J. Lauret and E. V. Shuryak, arXiv:nucl-th/0 110037. 

85. C. Nonaka and S. A. Bass, Phys. Rev. C 75, 014902 (2007) [arXiv:nucl-th/0607018]. 

86. D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 86, 4783 (2001) [arXiv:nucl- 
th/001 10581. 

87. D. Teaney, Phys. Rev. C 68, 034913 (2003) [arXiv:nucl-th/0301099]. 

88. P Romatschke and U. Romatschke, Phys. Rev. Lett. 99, 172301 (2007) [arXiv:0706.1522 
[nucl-th]]. 

89. M. Luzum and R Romatschke, Phys. Rev. C 78, 034915 (2008) [arXiv:0804.4015 [nucl-th]]. 

90. H. Song and U. W. Heinz, Phys. Lett. B 658, 279 (2008) [arXiv:0709.0742 [nucl-th]]. 

91. H. Song andU. W. Heinz, Phys. Rev. C 77, 064901 (2008) [arXiv:0712.3715 [nucl-th]]. 

92. K. Dusling and D. Teaney, Phys. Rev. C 77, 034905 (2008) [arXiv:07 10.5932 [nucl-th]]. 

93. H. Song and U. W. Heinz, Phys. Rev. C 78, 024902 (2008) [arXiv:0805.1756 [nucl-th]]. 

94. Z. Xu, C. Greiner and H. Stocker, Phys. Rev. Lett. 101, 082302 (2008) [arXiv:0711.0961 
[nucl-th]]. 

95. C. E. Aguiar, Y. Hama, T. Kodama and T. Osada, Nucl. Phys. A 698, 639 (2002) [arXiv:hep- 
ph/0106266]. 

96. M. MiUer and R. Snellmgs, arXiv:nucl-ex/0312008. 

97. S. Manly et al. [PHOBOS Collaboration], Nucl. Phys. A 774, 523 (2006) [arXivmucl- 
ex/0510031]. 

98. W. Broniowski, R Bozek and M. Rybczynski, Phys. Rev. C 76, 054905 (2007) 
[arXiv:0706.4266 |nucl-thl|. 

99. R. Andrade, F. Grassi, Y. Hania, T. Kodama and O. . J. Socolowski, Phys. Rev. Lett. 97, 
202302 (2006) [arXiv:nucl-th/0608067]. 

100. B. Alver et al. [PHOBOS Collaboration], Phys. Rev. Lett. 98, 242302 (2007) [arXivmucl- 
ex/0610037]. 

101. G. Wang, D. Keane, A. Tang and S. A. Voloshin, Phys. Rev. C 76, 024907 (2007) [arXivmucl- 
ex/0611001]. 

102. X. 1. Zhu, M. Bleicher and H. Stoecker, Phys. Rev. C 72, 064911 (2005) [arXivmucl- 

th/05090811. 

103. X. 1. Zhu, M. Bleicher and H. Stoecker, J. Phys. G 32, 2181 (2006) [arXivmucl-th/0601049]. 

104. B. Alver et al. [PHOBOS Collaboration], arXiv:0804.4297 [nucl-ex]. 

105. P Sorensen [STAR Collaboration], arXiv:0808.0356 [nucl-ex]. 

106. B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. C 72, 051901 (2005) [arXivmucl- 
ex/0407012]. 

107. R. SneUings [STAR Collaboration and ALICE Collaboration], Eur. Phys. J. C 49, 87 (2007) 
[arXiv:nucl-ex/0610010]. 

108. N. Borghini and J. Y. OUitrault, Phys. Lett. B 642, 227 (2006) [arXiv:nucl-th/0506045]. 

109. P Huovinen, Nucl. Phys. A 761, 296 (2005) [arXivmucl-th/0505036]. 

110. R. Snellings, arXivmucl-ex/03 10019. 

111. J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 93, 012301 (2004) [arXivmucl- 

ex/0312009]. 

112. D. Molnar and S. A. Voloshin, Phys. Rev. Lett. 91, 092301 (2003) [arXiv:nucl-th/0302014]. 

113. R. J. Fries, V. Greco and R Sorensen, arXiv:0807.4939 [nucl-th]. 



Collective phenomena in non-central nuclear collisions October 30, 2008 Draft 



55 



114. B. I. Abelev et al. [the STAR Collaboration], Phys. Rev. C 75, 054906 (2007) [arXivmucl- 
ex/0701010]. 

115. B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 99, 112301 (2007) [arXiv:nucl- 
ex/0703033]. 

116. A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 162301 (2007) [arXiv:nucl- 
ex/0608033]. 

117. M. Csanad et al, arXiv:nucl-th/05 12078. 

118. L. RavagU, H. van Hees and R. Rapp, arXiv:0806.2055 [hep-ph]. 

119. S. Pratt and S. Pal, Nucl. Phys. A 749, 268 (2005) [Phys. Rev. C 71, 014905 (2005)] 
[arXiv:nucl-th/0409038]. 

120. D. Mohiar and M. Gyulassy, Nucl. Phys. A 697, 495 (2002) [Erratum-ibid. A 703, 893 
(2002)] [arXiv:nucl-th/0104073]. 

121. R F. Kolb, P Huovinen, U. W. Heinz and H. Heiselberg, Phys. Lett. B 500, 232 (2001) 
[arXiv:hep-ph/00121371. 

122. R. J. M. Snellings, A. M. Poskanzer and S. A. Voloshin, arXiv:nucl-ex/9904003. 

123. M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 (1990). 

124. X. N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992). 

125. X. N. Wang, Phys. Rev. C 63, 054902 (2001) [arXiv:nucl-th/0009019]. 

126. M. Gyulassy, 1. Vitev, X. N. Wang and R Huovinen, Phys. Lett. B 526, 301 (2002) 
[arXiv:nucl-th/0109063]. 

127. M. Gyulassy, L Vitev and X. N. Wang, Phys. Rev. Lett. 86, 2537 (2001) [arXivmucl- 
th/0012092]. 

128. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 90, 032301 (2003) [arXivmucl- 

ex/0206006]. 

129. J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 93, 252301 (2004) [arXiv:nucl- 
ex/0407007]. 

130. Yuting Bai, PhD Thesis, NIKHEF, Amsterdam, Netherlands 

131. S. A. Voloshin [STAR Collaboration], J. Phys. G 34, S883 (2007) [arXiv:nucl-ex/0701038]. 

132. Z. w. Lin and D. Molnar, Phys. Rev. C 68, 044901 (2003) [arXiv:nucl-th/0304045]. 

133. S. S. Adler et al. [PHENK Collaboration], Phys. Rev. C 72, 024901 (2005) [arXivmucl- 
ex/0502009]. 

134. T. C. Awes [PHENK Collaboration], arXiv:0805.1636 [nucl-ex]. 

135. H. van Hees, V. Greco and R. Rapp, Phys. Rev. C 73, 034913 (2006) [arXivmucl- 
th/0508055]. 

136. D. Krieg and M. Bleicher, arXiv:0806.0736 [nucl-th]. 

137. R. Chatterjee, E. S. Frodermann, U. W. Heinz and D. K. Srivastava, Phys. Rev. Lett. 96, 
202302 (2006) [arXiv:nucl-th/051 1079]. 

138. W. Liu and R. J. Fries, arXiv:0805.3721 [nucl-th]. 

139. S. Turbide, C. Gale and R. J. Fries, Phys. Rev. Lett. 96, 032303 (2006) [arXiv:hep- 

ph/05082011. 

140. S. Huang [PHENIX Collaboration], arXiv:0804.4864 [nucl-ex]. 

141. C. Gombeaud and J. Y. OUitrault, Phys. Rev. C 77, 054904 (2008) [arXiv:nucl-th/0702075]. 

142. R R Kolb, L. W. Chen, V. Greco and C. M. Ko, Phys. Rev. C 69, 051901 (2004) [arXivmucl- 
th/0402049]. 

143. S. A. Voloshin and W. E. Cleland, Phys. Rev. C 53, 896 (1996) [arXiv:nucl- 
th/95090251;Phys. Rev. C 54, 3212 (1996) [arXiv:nucl-th/96060331. 

144. S. Voloshin, R. Lednicky, S. Panitkin and N. Xu, Phys. Rev. Lett. 79, 4766 (1997) 
[arXiv:nucl-th/9708044]. 

145. M. A. Lisa, S. Pratt, R. Soltz and U. Wiedemann, Ann. Rev. Nucl. Part. Sci. 55, 357 (2005) 
[arXiv:nucl-ex/0505014]. 

146. J. Bielcikova, S. Esumi, K. Filimonov, S. Voloshin and J. R Wurm, Phys. Rev. C 69, 021901 
(2004) [arXiv:nucl-ex/0311007]. 

147. H. Stocker, arXiv:0710.5089 [hep-ph]. 

148. C. NepaU, G. L Fai and D. Keane, Phys. Rev. C 76, 051902 (2007) [Erratum-ibid. C 76, 
069903 (2007)] [arXiv:0709.1497 [hep-ph]]. 



56 



S.A. Voloshin, A.M. Poskanzer, and R. Snellings 



149. N. Armesto etal, J. Phys. G 35, 054001 (2008) [arXiv:07 11.0974 [hep-ph]]. 

150. W. Busza, J. Phys. G 35, 044040 (2008) [arXiv:0710.2293 [nucl-ex]]. 

151. N. Borghini and U. A. Wiedemann, J. Phys. G 35, 023001 (2008) [arXiv:0707.0564 [hep- 
ph]]. 

152. D. Krieg and M. Bleicher, arXiv:0708.3015 [nucl-th]. 



